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FLAME  PROPAGATION 


CRITICAL  REVIEW  OF  EXISTING  THEORIES 


Ac  Introduction 

This  report  is  intended  as  a  comprehensive  review  of  the  present 
state  of  the  theory  of  flame  propagation  in  gases.  The  case  of  detona¬ 
tion  will  not  be  included,  the  discussion  being  limited  to  slow  combus¬ 
tion  (deflagration). 

The  main  purpose  of  the  theory  of  flame  propagation  is  to  relate 
the  velocity  of  the  flame  front  relative  to  the  unburnt  gas  to  fundamen¬ 
tal  physical  and  chemical  properties  of  the  combustible  mixture. 

Though  the  first  attempts  to  carry  out  this  proposition  were  made 
more  than  fifty  years  ago,  this  problem  is  still  far  from  being  solved, 
end  what  has  been  achieved  is  of  little  value  for  practical  application. 

The  main  reason  for  this  -failure  became  apparent  when  research  on 
the  kinetics  of  combustion  reactions  revealed  their  complicated  nature. 

It  is  now  common  knowledge  that  combustion  reactions  take  place  by  chain 
mechanisms,  involving  unstable  intermediate  products.  However,  it  has 
not  been  possible  yet  to  obtain  a  conclusive  knowledge  of  the  details  of 
combustion  reactions.  Without  this  knowledge  the  theory  of  flame  propa¬ 
gation  has  to  be  based  on  hypothetical  reaction  mechenisms;  this  fact  is 
primarily  responsible  for  the  discrepancies  between  theory  and  experiment. 

This  does  not  mean  that  the  theory  in  its  present  form  doe6  not  have 
considerable  value.  Apart  from  forming  a  basis  for  future  development,  it 


8 hove  that  the  importance  of  studying  flame  propagation  goes  beyond  the 
original  purpose  of  calculating  flame  speeds.  Once  a  definite  knowledge 
of  reaction  mechanisms  is  available,  the  theory  of  flame  propagation  will 
give  an  insight  into  the  complicated  interaction  of  chemical  reactions, 
diffusion  and  heat  conduction  taking  place  in  the  combustion  zone.  There¬ 
fore,  the  iaportance  of  further  research  in  the  field  of  kinetics  of  com¬ 
bustion  reactions  cannot  be  overemphasized. 

A  theory  based  on  the  effects  of  heat  conduction,  diffusion  and 
chemical  kinetics,  if  successful,  vould  enable  the  calculation  of  flame 
speeds  under  simple  flow  conditions,  e.g.  those  in  the  well-known  soap- 
bubble  and  Bunsen  burner  methods.  A  more  complete  theory  would  have  to 
include  also  the  modifications  introduced  by  more  complicated  flow  condi¬ 
tions,  particularly  by  turbulence.  The  attempts  to  include  these  effects 
in  the  theory  are  only  at  their  beginning;  the  experimental  evaluation 
of  this  aspect  of  flame  propagation  likewise  presents  difficulties  not 
yet  overcome. 

Another  complication  arises  in  many  technical  applications  because 
the  assumption  of  a  homogeneous  mixture  of  fuel  and  air,  on  which  the 
usual  theories  of  flame  propagation  are  based,  is  not  fulfilled.  A 
special  case,  namely,  combustion  in  the  mixing  zone  of  fuel  and  air,  with 
mixing  taking  place  exclusively  by  diffusion,  has  been  treated  theoreti¬ 
cally  with  remarkable  success.  Future  development  will  have  to  include 
the  cases  of  convective  mixing  and  of  partially  premixed  gases  and  should 
establish  thus  the  connection  with  the  theory  of  propagation  in  homogeneous 
mixtures. 
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Bo  Definition  of  Transformation  Velocity 


Before  discussion  the  theory  of  flame  propagation,  it  is  necessary 
to  find  a  definition  of  the  velocity  of  flame  motion  vhich  is  unique  and 
independent  of  the  motion  and  local  state  of  the  v  '  ^rnt  gas. 

One  commonly  encountered  definition  la: 

"Flame  velocity  is  the  velocity  of  the  flame  front  in  a  direction 
perpendicular  to  its  surface,  relative  to  the  unburnt  gas  where  it  is  at 
initial  conditions." 


A  more  satisfactory  definition  is  the  following: 

"The  flame  velocity  is  equal  to  the  mass  flow  of  unburnt  gas  through 
unit  area  of  the  flame  front,  divided  by  the  initial  density." 


Under  stationary  conditions  of  one-dimensional  flow,  the  second 
definition  is  equivalent  to  the  first  one. 


This  can  be  demonstrated  by  the  relations  in  the  Bunsen  burner 

cone,  which  is  shown  in  the  idealised  diagram  Figure  1.  V  is  the  flow 

velocity  of  the  unburnt  gas,  Vn  it's  component  perpendicular  to  the  cone 

surface,  dOlan  element  of  cone  surface  and  do  its  projection  on  the  base. 

The  mass  flow  through  dUJ  is  given  by  dqs^vdOs  ^  V  dui  sin  OC 

where  J*9  is  the  initial  density.  Because  of  the  relation  Vn  *  V  sin  OC 

there  follows  Vn  =  — 

n  J>9  du) 


As  Vp  corresponds  to  the  first  definition,  and  tj— 
definition  of  propagation  velocity,  the  equivalence  of  the  two  is  thus 
proven. 


dq 

dU> 


to  the  second 
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It  is  impossible  to  observe  this  velocity  directly,  as  conditions 
where  the  unburnt  gas  is  at  rest  cannot  be  realized  experimentally.  The 
velocity  of  propagation  must  be  calculated  either  as  the  vector  differ¬ 
ence  of  the  observed  velocity  of  flame  movement  and  velocity  of  flow  of 
the  unburot  gas,  or  more  convenitntly  as  total  mass  flow  through  the 
flame  front,  divided  by  the  surface  of  the  flame  front  and  by  the 
initial  density , 

It  should  be  mentioned  that  unfortunately  a  confusing  variety  of 
terns  is  used  in  the  literature  on  flame  propagation.  In  this  report 
the  term  "transformation  velocity"  ,  S ^ ,  will  be  adopted  for  the  quantity 
defined  above,  in  accordance  with  the  National  Bureau  of  Standards  and 
the  NACA1. 


Co  The  Thermal  Theories 

(1)  The  Older  Thermal  Theories 

The  thermal,  theories  are  based  on  the  equation  of  balance  of 
thermal  energy  in  the  combustion  zone.  The  objections  raised 
by  various  autho^8^,^,^,^,®,<^,^a,^,^  against  the  purely 
thermal  point  of  view  in  treating  flame  propagation  will  be 
discussed  in  later  chapters. 


*  Other  terms  frequently  used  for  the  same  quantity  are:  burning 

velocity2,  fundamental  flame  speed*3,  flame  speed*,  ignition  velocity^, 
and  normal  velocity^8. 
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The  equation  of  heat  balance  can  be  derived  es  follows  :  de¬ 
noting  by  H  the  heat  content  of  unit  mass,  by  J>  the  density 
of  the  gas,  by  jfl  the  vector  of  heat  flow,  by  W  the  reaction 
rate  and  by  Q  the  heat  of  reaction,  referred  to  unit  mass  of 
combustible,  the  rate  of  change  of  heat  content  per  unit  volume 
is  given  by 


=  -dl»7f+  wQ 

The  heat  flow  vector  is  given  by 


(1) 


£1  s  -  X  grad  T  +  j)  Hv  (2) 

where  X  is  the  heat  conduction  coefficient,  T  the  absolute 
temperature,  and  V  the  velocity  of  gas  motion.  Combining 
equations  (1)  and  (2)  with  the  continuity  equation 


1L, 


di*(/v) 


there  results  (taking  in  account  that 


(3) 


dlv(/Hv)  s  H  dlv  (  /V)  +  J>V  grod  H) 

div(Xgrod  T)  —  jPVgrod  H  +  wQ  (4) 

As  the  treatment  is  concerned  with  flame  propagation  in  uniform 
medium  under  stationary  conditions,  the  coordinate  system  can 
be  fixed  in  the  flame  front.  This  causes  all  time  derivatives 
to  vsnish.  By  disregarding  the  curvature  of  the  flame  front, 
the  equations  can  be  specialized  further  to  the  case  of  one¬ 
dimensional  flow.  The  continuity  equation  under  these  condi¬ 


tions  becomes 


)  =  0 
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where  X  is  the  direction  of  flow,  which 


leads  to 


jdV  !  q  *  f>VQz  Con»t.  (5) 

where  q  ia  the  mass  flow,  the  initial  density,  and 
V0  the  velocity  of  the  unburnt  gee  at  initial  conditione 
relative  to  the  flame  front,  which  by  definition  is  identical 
with  the  transformation  velocity  .  Assuming  in  addition  X 

Jii  JT  _ 

to  be  independent  of  temperature,  and  replacing  -SE  by  C0~* 

dx  ^  dx 

where  Cp  is  the  specific  heat  at  constant  pressure,  equation 
(4)  becomes  finally 


d*T  AC  P 
dX*  X 


(6) 


The  problem  consists  in  finding  that  value  of  which  is 
compatible  with  a  solution  T{X)  of  equation  (6)  which  fulfills 
the  boundary  conditions: 

T  =  T0  for  X  •  —  oo 
T  =  Tf  for  X  =  + 

where  T#  and  Tf  are  the  initial  and  final  temperatures  re¬ 
spectively  „ 


*  This  approximation  can  be  shown  to  be  reasonably  valid  (ref.  6,  footnote 
p.  20) $  with  the  same  degree  of  approximation,  Q  can  be  assumed  to  be 
independent  of  temperature. 
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The  main  complication  lies  in  the  dependence  of  the  reaction 
rate  Uf  on  the  temperature  and  concentrations  of,  in  general, 
a  number  of  aubatoncee  participating  in  the  reaction.  These 
concentrations,  as  functions  of  x  are  determined  by  diffu¬ 
sion,  generation  and  destruction  of  particles  according  to 
the  particular  reaction  mechanism.  The  lack  of  knowledge  of 
the  reaction  mechenisma  of  the  most  important  combustion  pro¬ 
cesses  precludes  at  present  an  exact  treatment  of  the  problem, 
as  will  be  discussed  more  in  detail  in  Chapter  G.  Only  by 
Introducing  rather  radically  simplifying  assumptions  can  an 
approximate  solution  be  found. 

The  older  attempts  of  solution  are  based  on  the  assumptions, 
that: 

(a)  Up  to  a  cerU-in  ignition  temperature,  Ti  which  is  taken 
to  be  a  unique  function  of  composition  end  pressure,  no 
reaction  tekes  place, 

(b)  Once  this  temperature  has  been  reached  the  reaction  pro¬ 
ceeds  at  a  constant  rate,  determined  only  by  the  initial 
conditions.  Assumptions  of  this  nature  led  to  the  follow¬ 
ing  expressions  for  the  transformation  velocity 

Tf-Tl  in 

Sf  =  K  ~=~  - - -r-  (Mallard  and  Le  Chetelier  u) 

1  Cp  Tj  -  T0 

S*  z  .  /  — — - i-  (Jouget,  Crussard*1) 

V-P0Cp  Tj  Tj  -T0 
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(Denlell12) 


St  = 
St  . 


/XJ  T„ 
V  A  cp  tn 


/ 


K  Xltf 

D  5p  0o 


Tf-T0 

T)-T. 

tFTT 

Tj  -  T0 


(Damkohler*^) 


tore  K  denotes  a  constant,  f  a  quantity  proportional 
to  the  reaction  rate,  Cp  the  average  specific  heat 

at  constant  pressure,  X  the  heat  conductivity,./^  and 
To  the  initial  density  and  absolute  temperature,  Tj 
the  ignition  temperature,  Tf  the  final  temperature 
(absolute)  and  a  harmonic  mean  teaperature  nearer  to 
Tf  than  to  TjUfthe  mean  reaction  rate  and  a0  the  initial 
concentration  of  combustible.  By  further  specializing  the 
formulas  for  H2  -  02  and  H2  -  0?  -  N£or  H?-Q,-CC^  jalxtures, 
the  following  expressions  were  derived-. 


Ap„  T02  Tf  -  Tj 
cp  Tj  -T. 


[^2]  [^2]  (Nuacelt^) 


SlS  y/K-H/il-H,)* 

j  \ - 


~y. 


"reaction-kinetic1'  term  "thermic  effect"  term  "initial  condition"  te 
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where  and  Jbgj  are  moler 
concentrations,  H2  the  raol  fraction  of  tydrogen,  OC  the 
mol  fraction  of  oxygen  in  the  atmosphere  Oj+N^or  Og+CO^h! 
the  heat  of  combustion  of  unit  volume  of  tydrogen,  Po 
Initial  total  pressure.  It  will  now  be  outlined  briefly 
how  an  expression  similar  to  those  cited  is  derived^,, 

The  origin  of  the  coordinate  system  is  placed  at  the  bound¬ 
ary  between  boating  and  reaction  zone,  the  positive  X  -axis 
pointing  towards  the  reaction  zone  (Figure  2) 

Thus,  for  X^O,  equation  (6)  reduces  to 


d*T 

dX* 


-OC  -~“s°  00  s  ^-*St(7) 

dx  x  x 


With  the  appropriate  boundary  conditions,  the  solution  is 


oc  X 


T-T0s(Tj-T#)-e  for  -  oo  X  0  (8) 


In  the  reaction  zone  of  thickness  dr  the  reaction  rate  W 
is  assumed  to  be  constant.  Tnus,  there  can  be  written 


W 


da  _  do  dx.Oo 
dt  dx  dt  "  dr 


St  (9) 


where  0  is  the  concentration  of  combustible,  with  the 
initial  value  Oo  Furthermore,  the  heat  of  reaction  can 
be  expressed  through  the  total  tempera turo  rise: 


o  _ 


(10) 


Q  ~p~  =  L  =  Cp  (Tf-T.) 

-'o 

where  L  is  the  heat  of  reaction  referred  to  unit  mass  of 
the  whole  mixture*  Introducing  sxpre6sions  (9)  and  (10),  the 
solution  of  the  fundamental  equation 


H2T 

-2-1  -  OC  +-=^3-  =0 

dx*  dx  X 


(6) 


for  the  reaction  zone,  O^X  ^dr  becomes 


T-t«(Tf-T.)Jf-  +(TrT0) 


l  -e 


ocdj- 


(ID 


In  order  to  determine  Sj  (which  is  contained  in  OC  ),  an 

additional  boundary  condition  must  be  introduced,  namely,  the 

continuity  of  ^1  at  the  boundary  of  heating  and  re-ction 

zone*  Equating  )x  from  equation  (8)  to 

(^X  \  from  equation  (11),  the  transcendental  equation 

dX  /X  =0 

for  OC 


Tj  -T0 
Tf-T. 


-ocd  r 

1  -  e 


oc  dr 


(12) 


is  obtained* 


10  - 


X  =  0  X  —  dr 

FIG.  2 

TEMPERATURE  DISTRIBUTION  IN  THE  FLAME  FRONT 


In  the  older  theories  the  ignition  temperature  Tj  was  identified 
with  the  one  determined  by  "static"  ignition  tests,  in  which  it 

was  found  to  be  much  lower  than  T*  .  Therefore.  ZLzZsL  was  assumed 

»  T  -  T 

f  *° 

to  be  small,  or  OCdr  large 


This  led  to  the  approximation 


Tf-T„ 

Wo 


-S  CP 

x 


•dr  •  St 


(13) 


Substituting  dr  from  equation  (9)  leads  to  the  final  result 


S  -  X  *  Tf  ~  To  (U) 

V  CP°°  Ti  T° 

This  is  ecsenti6lly(  Identical  w.ith  the  expressions  given  on  pages  7  &  8 
The  various  forms  of  the  expressii  r.s  result  from  special  assump¬ 
tions  about  the  dependence  of  the  reaction  rate  on  the  initial 
conditions . 

(2)  The  Theory  of  Zeldovich  and  Frank-KamenetBky^,^,^'^ 


Tne  older  thermal  theories  have  been  criticized  by  various 
authors  on  the  ground  that  they  are  based  on  assumptions 
which  cannot  be  justified  by  the  current  theory  of  reaction 
kinetics  Theoretical  considerations  and  a  great  amount  of 
experimental  evidence 9,16,17  established  that  the 
temperature  of  ignition  of  an  explosive  gaseous  mixture 
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depends  on  other  variables  beside  initial  composition  and  pres¬ 
sure  <>  Most  important  for  the  theory  of  flame  propagation  is  the 
relation  between  the  ignition  temperature  and  the  time  lag 
(induction  period)  of  ignition.  Whereas  in  the  usual  experi¬ 
mental  determinations  of  ignition  temperatures  this  time  lag  is 
of  the  order  of  0.1  to  several  seconds,  in  a  flame  travelling  at 

only  1  m/sec.  and  having  a  thickness  of  the  heating  zone  of  0.1 

* 

mm  the  induction  period  is  only  10  ^  sec.  From  this  considera- 

6>  V  18  IS 

tion  the  authors  of  the  newer  thermal  theory  *  *  *  conclude 


that  the  ignition  temperature  in  a  flame  must  be  taken  very  high, 

approaching  the  final  temperature  .  More  explicitly,  it  is 
7  25 

stated  '  ,  that  it  can  be  lower  than  only  by  an  amount  of  the 

order  of  -Tf-,  where  E  is  the  activation  energy  of  the  reaction 

E 


and  R  the  gas  constant. 


This  conclusion  is  in  agreement  with 


the  experimental  fact,  that  the  transformation  velocity  does  not 


reach  infinity  on  raising  the  initial  temperature  to  values  as 
high  a6,  or  even  higher  than,  the  "static"  ignition  temperature, 
as  would  follow  from  the  expressions  derived  according  to  the 


older  theories.  In  accordance  with  these  considerations  which 


show  that  the  ignition  temperature  Tj  is  not  a  well-defined 
physical  quantity,  Zeldovich  and  Frank-Kamenetsky  introduce  it 
merely  for  convenience  of  computation;  no  particular  significance 
is  attributed  to  it  and  it  is  eliminated  from  the  final  expressions. 
The  reaction  rate  is  introduced  now  as  a  function  of  temperature 
and  concentration.  The  derivation  as  given  by  Semenov  will  be 


followed  here. 
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As  a  first  step,  an  approximate  solution  of  equation  (6)  for 


the  reaction  zone,(o^X^<i ,  Tj^T^Tf)!®  derived  by  neglecting 
the  terra  g  ^  — L.  .  (Semenov  justifies  this  uimplifi- 

cation  by  taking  into  account  equations  (9)  and  (10) j  since 
equation  (9)  is  based  on  the  assumption  of  a  constent  reaction 
rate,  Semenov's  proof  cannot  be  considered  correct.)  Neglect- 

and  in- 


ing  the  term  with  -  T~  replacing  — by  -il 

dx  dx  dx 


d 

dT 


tegrating  equation  (6)  over  the  reaction  zone  (o  ^X  ^dr  ) 
there  results 


Ub 


■^-1-  vanishes  for 

dx 


x  -  dr) 


(15) 


For  —  oo^X  ,  solution  (8)  is  carried  over  into  the  new 
theory.  This  gives 


('5r)xr0oc  (t  -TJs°c(Tf-T*>  <16> 

The  last  approximation  is  justified  since  Tj  is  now  assumed  to 
be  close  to  T|  .  Substituting  for  (T^— T0)from  equation  (10), 
and  equating  (15)  with  (16),  the  general  expression  for  the 
transformation  velocity  is  obtained; 


-  13  - 


A  suitable  lav  for  the  reaction  rate  oust  now  be  Introduced 
Into  this  expression.  Here  the  main  difficulty  arises,  due  to 
the  complex  character  of  practically  all  the  important  combus¬ 
tion  reactions.  Zeldovich  and  Frank-Kamenetsky  have  applied 
the  new  theory  only  vhere  a  "classical "  (Arrhenius)  kinetic  lav 
can  be  assumed  •  (For  comprehensive  discussions  of  chemical 
kinetics  in  gases,  see  refs.  20,  21) 

Thus,  there  is  introduced: 


w-’S-e 


w  =  Ka  e 


E 

“rT 


W  =  Ka2e 


for  a  zero  order  reaction 

for  a  1st  order  reaction  (18) 

for  a  2nd  order  reaction 


Here,  E  is  the  activation  energy,  a  the  concentration  of  com¬ 
bustible,  and  the  factors  S  and  K  depend  essentially  only  on 

the  collision  frequency  of  the  molecules,  which  can  be  calculated 

20  21 

by  the  methods  of  statistical  mechanics  *  . 


-#r 


The  factor  ©  increases  very  rapidly  with  temperature  for 

the  usual  values  of  the  activation  energy  (25,000  to  80,000  cai/mol) 


*  This  does  not  mean  that  the  mechanism  of  the  reaction  is  assumed  to  be 
simple  and,  respectively,  of  zero  order,  aono-or  biaolecular.  Compare 
ref.  21  for  this  essential  distinction  between  order  and  mechanism  of 
a  reaction. 
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Therefore,  the  temperature  dependence  of  S  or  K  is  neglected, 
their  value  being  taken  at  the  highest  temperature  Tf  .  Since 
the  difference  0=Tf-T  ia  small  compared  with  T,  the  follow¬ 
ing  approximation  can  be  made: 


-  6 


Rv 

e  t 


In  the  case  of  sero  order  reaction,  the  integral  I  in  equa¬ 


tion  (17)  can  now  be  evaluated  immediately: 

fT »  -if,  ' 6  "If.  RTf*  rP I 

X=  /  Se  RT<  •  •  RT»  dT -  Se  Jt- 

Ti  r  .  ^  0 


with 


It  is  seen  that  the  integral  r 

converges  towards  the  value  one  very  rapidly  for  increasing^; 
The  authors  of  this  theory  set  the  integral  equal  to  one,  on 


the  assumption  that  the  term 
justify  both  conditions: 


is  always  large  enough  to 


j  CL  I  equivalent  to  >  2  ,  »nd 

Tf  -Tj  ^  S  Tf  ,  which  has  been  previously  assumed  to  be 
valid  (p0  12),  The  remarkable  result  of  this  assumption 

is  the  elimination  of  the  ignition  temperature  from  the  final 
expression. 
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The  formula  for  the  transformation  velocity  for  aero  order 
reaction  becomes  thus 


S  ->  /2.^-  e,~  £ff  £2*  _v  Za  i* _ Tfff  RTf‘ 

St=Vj»o.LS*  *  E  “Wp  0.  E(T,-T<J 

A  new  problem  arises  in  the  treatment  of  first-  and  second- 
order  reactions |  because  the  concentration  0  has  to  be 
brought  into  a  functional  relation  to  the  temperature.  By 
transforming  equation  (6)  suitably,  and  introducing  the 
diffusion  equation  for  a  it  can  be  "  iown^ 
iables 


that  the  var- 


and  6  s 


(T-TJ 


obey  the  differential  equations 

x  dae  ae 


ep  dx* 


'  S*dx  +  ws° 

of&r- st$r  +  "  =  ° 

where  D  is  the  coefficient  of  diffusion 


(21) 


Both  variables  have,  in  addition,  to  fulfill  the  same  bound¬ 
ary  conditions: 


(o'O-it0).  io0 


(a*; 


(«Mjr  =  -5-  =  ^(T rT.Me)+, 

If  the  following  condition  is  satisfied: 


0  =  JTP=I  (®) 

where  X  i8  the  coefficient  of  temperature  conduction,  OC 
must  be  Identical  with  G  throughout  the  whole  range 
<*<+ 


(20) 
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Condition  (22)  is  never  exactly  fulfilled;  however,  the  elemen¬ 
tary  kinetic  theory  of  gosos  gives  as  a  first  approximation 
Xs  y>cv  and  a  better  approximation  leads  toX  =  K^>CVD, 
with  Kcy  —  Cp  (ref.  22) 


This  result  can  be  written 


cpT  +  *  CpT.  +  -%T  «  CpTf  =  Con,t. 

aO 


(23) 


Since  c.T  is  the  enthalpy  and  u  —  the  chemical  energy  avail- 

P  / 

able  per  unit  volume,  (23)  expresses  the  constancy  of  the  sum 
of  thermal  and  chemical  energies  during  the  whole  coabustion 
process c.  This  relation  had  already  been  postulated  without 
proof  by  Lewis  and  von  Elbe2,  in  their  theory,  which  will  be 
discussed  later.  It  is  very  important  to  note  that  it  holds 
only  if  condition  (22)  is  fulfilled. 

Equation  (23)  can  be  transformed  by  meanB  of  (10)  to 


a  =  a. 


Jill  SO  Jt 

T.  -T.  °«  T 


T,-T 


TfT. 


(23' ) 


4  '( 

On  evaluating  the  integral  I  in  (17),  the  approximations  (19) 
are  again  introduced,  and  T  in  the  denominator  of  (23')  is 
approximately  put  equal  to  Tf.  Tj  is  eliminated  again  in  an 
analogous  way  ao  in  the  case  of  zero  order  reactions.  The  ex¬ 
pressions  for  the  transformation  velocity  become  thus: 
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(T)‘ 


E 

•*T, 


2XcpK 


for  first  order  reaction  ( W  =  K  0  6 


) 


and 


(24) 


E 


for  second  order  reaction  (w  =  Ka2e  ^  ) 

IQ 

In  one  of  their  papers  '  Zeldovich  and  Frank-Kamenetsky 
give  a  general  expression  for  reactions  of  order  ps 


St 


P-M 


(26) 


where  X  is  the  temperature  conductivity  ( ~  — “—JjTfthe  tine 

>’CP 

during  which  the  reaction  would  proceed  isothermally  at  the 
temperature  Tf  and  at  the  initial  concentration  0o 


( 


Tf 


=  Ka 


P-i 

0 


and 
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0f  is  equal  to  - \  (Tf  —  T0).  Equations  (20),  (2/*)  and 

RT<  T 

(2 5)  are  seen  to  ogre®  with  1 26) ,  except  for  the  factors 

/T\*  T« 

I — “-),  respectively,  which  enter  the  expressions  due  to  a 

V  Tf  / 


ond 


correction  on  the  diffusion  equation  applied  by  Semenov 
(See  ref.  6,  foot  notes  pp.  35  and  41). 

It  is  evident  that  the  effects  of  temperature  and  concentration 
on  the  reaction  rate  on  cessing  through  the  combustion  zone 
oppose  each  other.  Semenov^  therefore  investigates  the  validity 
of  the  assumption  previously  made  (see  p  12  that  practically  no 
reaction  takes  place  until  the  temperature  approaches  the  final 
temperature  Tf  .  For  the  case  of  second-order  reaction,  he  shows 
that  the  assumption  Tf  -  Tj  =  0.25Tf  is  valid  only  for  values  of 
^  S  0.1  or,  with  Tf  =  2000°,  for  E  £=  40,000  cal/mci* 


Obviously,  for  the  case  of  first-order  reaction,  a  less  rigorous 
restriction  applies. 


Formulas  (24)  and  (25)  can  be  formally  obtained  by  introducing 
into  the  expression  for  zero-urder  reaction,  (20),  an  effective 
concentration 


a  eft  = 

and  setting  S=  Ko 


°.CP  /RT,*\  Jo  J> Cp  RTf  T„ 

L  \  E  /  Tf  '  0  E  Tf 

2 

cfffor  first-order  andS=KOgff 


(27) 

for  second-order 


*  Hb  points  out,  in  this  connection,  that  a  factor  of  2  instead  of  4  under 
the  square  root  in  equation  (25)  is  more  nearly  correct. 
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reactions.  This  allows  the  extension  of  the  theory  to  cases 
where  the  reaction  rate  depends  on  the  concentrations  of  more 
than  one  chemical  species.  For  instance,  if  the  kinetic  law 

-Jr 

W  sK  a  b't  n  can  be  applied,  and  there  is  a  deficiency  of  a 
and  an  excess  of  b  ,  b^ff  may  be  set  equal  to  its  concentra¬ 
tion  in  the  combustion  products  and  calculated  according 

to  (27) o  Semenov^  calls  attention  to  the  second  form  of 
equation  (27),  which  shows  that  Oeff  does  not  depend  directly 
on  the  initial  concentration  0o  ;  this  circumstance  imposes 
considerable  limitations  on  the  possibility  to  draw  conclusions 
on  reaction  kinetics  from  flame  propagation  data. 

Finally,  the  theory  is  extended  to  cases  where:  (1)  specific 
head  and  heat  conductivity  vary  with  temperature;  (2)  the 
nimber  of  molecules  changes  during  reaction,  and  (3)  the  diffu¬ 
sion  coefficient  is  not  equal  to  the  coefficient  of  temperature 
conduction  ^condition  (22)J„  The  final  expressions  given  are: 


for  first-order  reaction,  and 


(24a) 


(25a) 
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Where  the  values  of  have  to  be  taken  for  the  reaction 
products  at  the  final  tempers t'ire  Tf  ,  D  is  the  diffusion  co¬ 
efficient  of  the  reacting  substance  in  the  reaction  products 
at  the  temperature  Tf  ,  hj  is  the  number  of  molecules  before 
and  n£  the  number  of  molecules  after  reaction. 


D „  Comparison  of  the  thermal  theories  with  experimental  results 
A  large  amount  of  experimental  material  on  flame  propagation 
velocities  exists,  and  many  attempts  have  been  made  to  verify  the  pre¬ 
dictions  of  the  older  thermal  theories  with  these  data.  It  should  be  re¬ 
membered  in  this  connection  that  the  older  tneories  do  not  allow  the  cal¬ 
culation  of  absolute  values  of  the  transformation  velocity,  as  the  ex¬ 
pressions  contain  an  unknown  constant  (  K  or  f  in  the  formulas  on  pages  7  &  3) 
For  some  particular  Initial  condition  in  a  series  of  experiments  the  theo¬ 
retical  value  can  thus  be  fitted  to  the  experimental  result..  The  incorrect 
assumptions  of  the  older  theories,  namely,  co'rtant  ignition  temperature 
and  constant  rate  of  heat  liberation  throughout  the  combustion  process,  are 
largely  responsible  for  the  numerous  discrepancies  between  experimental  and 
theoretical  values.. 


Formulas  of  the  form 


TlzI 


■t  w  vrf:  (r,f-l3) 

allow  the  prediction,  in  a  qualitative  way,  of  the  existence  of  limits  of 
inflammability,  as  for  very  lean  or  very  rich  mixtures  Tf  must  become  equal 
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to  Tj  and  must  become  zero.  Qualitative  predictions  of  the  depend¬ 
ence  of  on  the  composition  of  the  gas  mixture  in  fairly  good  agreement 
with  experimental  results  were  made  by  Jahn^u 

IQ 

Lewis  and  von  Elbe*7  have  recently  emphasized  that  the  formula  of 
Mallard  and  Le  C  ha  teller^  in  the  form 

s.  =  Tf~T* — !— 

1  J>  Cp  Tj-To  dr 

offers  possibilities  for  the  study  of  combustion  which  have  not  been  fully 
appreciated  in  the  past.  They  propose  a  semi-empirical  approach,  combin¬ 
ing  experimental  values  of  Sf  and  dr  with  the  Mallard-LeCha teller 
treatment  for  deriving  other  useful  information,  such  as  the  overall  re¬ 
action  rate.  They  assert  also  that  the  rather  difficult  measurement  of 
dr  can  be  replaced  by  the  determination  of  the  more  easily  measurable 
quenching  distance, that  is  the  distance  between  the  flame  front  and  a  solid 
heat-conducting  surface^  A  theory  linking  the  minimum  spark  energy  for 
ignition  with  quenching  distance  and  transformation  velocity  was  also  pro¬ 
posed  by  these  authors.  For  a  more  detailed  critical  discussion  of  the 

oldor  thermal  theories,  reference  is  made  to  the  reviews  by  Lewis  and  von 

2  3  8  9 

Elbe  ,  Coward  and  Payman^  and  Jost  *  . 

The  new  thermal  theory  has  been  applied  so  far  only  to  the  com¬ 
bustion  of  carbon  monoxide,  and  of  the  liquid  nitroglycol^.  This  restric¬ 
tion  is  due  to  the  fact  that  the  authors  of  this  theory  doubt  its  appl  ca- 
bility  to  fast  reactions  of  the  branched  chain  type,  particularly  to  the 
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combustion  of  hydrogen  and  the  hydrocarbons  Though  the  actual  re¬ 
action  mechanism  of  the  carbon  monoxide  oxidation  is  of  a  chain  character 
and  its  details  are  not  established  conclusively,  a  kinetic  law  of  the 
general  form  (18), 


w  .  -  .  -2^1  .  Igs]  .  Ze-«T  •  [co]  VJ 


(18') 


with  E  =  25000  col,  is  assumed  by  Semenov  for  this  reaction. 

This  law  takes  care  of  the  observation  that  moisture  has  a  profound  in¬ 
fluence  on  the  flame  propagation  in  CO  (cf.ref.l)  Two  cases  have  to 
be  distinguished,  on  introducing  (18')  into  equation  (17): 


(1) 

CO 


Excess  of  CO  in  the  initial  mixture.  Here  the  concentrations  of 
and  HgO  are  both  taken  at  the  final  temperature  and  in  the  burnt 


gas,  by  multiplying  each  with  -^2.  and  setting  M,  •  [co]0  ■  p°Jf 


Thus  (18:)  becomes  W  =  2  e 


RT'  {%)'  '  [Ht0]o  [C0]f 


The  transformation  velocity  is  computed  by  applying  equation  (17), 
with  the  following  substitutions: 


Zs  hT  ’  vt^cp  (2: 

where  Jl  is  the  mean  free  path,  C  the  mean  molecular  velocity 
and  n  the  number  of  molecules  per  unit  volume  of  the  gas  at  the 
temperature  Tj  ■  By  furthermore  setting 


_Q_  = 


LPo  py, 


0  [codf  ' 
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1b  the 


where  ie  the  molecular  weight  of  CO  and  0o  =  n  -^r 
number  of  molecules  per  unit  volume  at  the  initial  temperature  TOJ 
and  by  assuming  the  numerical  values 

Q  =  7.I04  “‘/mol,  a-  30  ,  cJ00  =  5.I04  cm/i«c 
the  following  expression  for  the  transformation  velocity  is  obtained: 


St=  /|400  -t— -® —  «  (?9) 

v  E  « 

(2)  Excess  of  oxygen  According  to  what  was  said  previously,  equation 
(27)  must  be  applied  in  this  case  to  introduce  the  effective  concen¬ 
tration  of  CO.  To  give  the  final  formula, 

[co],„  .  SBA  llL 

has  to  be  substituted  for  fco]^  in  (29) 

Without  giving  numerical  results,  Semenov^  states  that,  the  agreement  of 
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these  expressions  with  the  experimental  date  of  Fiock  ,  Jahn  ,  Passauer 
and  unpublished  observations  by  Ear  sky  is  satisfactory  Figure  3  shows  a 
plot  taken  from  Semenov’s  paper  showing  a  comparison  of  his  results  with 
the  data  of  Chitrin,  Barcki,  and  Jahn*  In  particular,  Barsky's  experi¬ 
ments  are  stated  to  confirm  the  following  predictions  of  the  theory: 

(a)  for  excess  oxygen  the  transformation  velocity  does  not  depend  on  either 
the  initial  CO  or  02  concentration,  if  by  suitable  inert  admixtures  the 
gas  composition  is  adjusted  in  such  e  way  that  the  final  temperature  Tf 
is  held  constant;  (b)  for  excess  CO,  the  transformation  velocity  varies 


with  /[C0]f  ^  again  Tf  is  maintained  constant.  It  is  further  stated 


0  20  40  60  %  CO 

FIG.  3 

TRANSFORMATION  VELOCITY  OF  CARBON  MONOXIDE- AIR  MIXTURES 


(AFTER  SEMENOV  REF.  6  ) 

FULL  CURVE  !  DATA  OF  PASSAUER 
X  !  DATA  OF  CHITRIN 
0  !  DATA  OF  BARSKY 
DASHED  CURVE  !  THEORETICAL  VALUES 
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s, 

CM/SEC 

600 


200 


0  I  2  Vh70% 

FIG.  4 

SPATIAL  FLAME  VELOCITY  VS.  SQUARE-ROOT  OF 
H,0  CONCENTRATION 

(REPLOT  OF  DATA  BY  FIOCK  AND  MARVIN 
REF.  I  BY  SEMENOV,  REF  6  ) 

1  :  Fc0  f  o  *  100  MM 

2  •  PcO  +  o[  z  ‘50  MM 

3  •  PC0  +  02  s  200MM 

4  !  Pco  +  0£  *  300 MM 
5!  PcO  +  o£  s  760MM 


i 


that.  Flock' r.  data  (ref.  1,  Fig.  9)  give  approximately  straight  if 

replotted  &g  S-*  v*.  ,  Figure  4  in  agreement  with  the  theo: 


replotted  &g  S#  v*.  v/[hso]  ,  Figure  4  in  agreement  with  the  theory 

Finally,  in  the  care  of  CO  -air  mixtures,  the  theory  predicts  independence 
of  the  transformation  velocity  from  the  total  pressure  for  constant  mois¬ 
ture  percentage  end  inverse  proportionality  to  v/P  for  constant  partial 
pressure  of  water.  This  is  said  to  be  in  agreement  with  experiment.  For 
C0-02  mixtures,  temperatures  around  3000°  K  are  reached  in  the  fit  me 
and  therefore  dissociation  must  be  taken  into  account.  Semenov  indicates 

that  the  theory  then  predicts  an  increase  of  with  p0**, while  experi- 

_0.£ 

ments  showed  an  increase  with  r 


The  application  of  the  thermal  theory  to  the  combustion  of  nitro- 
flycol  will  not  be  discussed  here,  as  combustion  on  the  surface  of  a  liquid 
is  beyond  the  scope  of  this  report.  A  similar  theory  of  flejne  propagation 
in  a  colloidal  propellant"^  should  also  be  mentioned  here,  without  further 

discussion. 


E .  Thickness  of  the  flame  front and  shape  of  flame s  near  walls  according 
to  the  newer  thermal  theory. 

Once  the  transformation  velocity  is  determined,  it  is  possible  to  com¬ 
pute  the  thicknesses  of  the  heating  zone  and  reaction  zone.  From  equations 
(v)  and  (3),  it  is  seen  that  the  temperature  difference  T~T0  drops  from 


T|-T. 


(Tj-To)  In  a  distance 


CP J> 

o 


♦Spatiel  flame  velocity  in  soap-bubble  experiments 
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This  vblue  i£  taken  as  a  measure  of  the  heating  zone  thicknese.  From 


equation  (16),  by  putting  the  mean  value  of 

and  by  integrating,  one  obtains 


d I 

dx 


equ»l  to  ±  (iL)s 


Ta  -T; 


r  =2-J - L  dt  ~  0.4  d* 

Tf-Te  1  1 


as  thicknees  of  the  reaction  zone.  With  the  usual  values  of  transforma¬ 


tion  velocity  St  and  temperature  conductivity  X,  d^  is  of  the  order  of 
0.2  to  0„3  nun,  end  dr  of  the  order  of  0.1  mm.  Semenov  deducee  further 

-5 

that  the  gas  molecules  stay  about  3 >10  sec  in  the  reaction  zone  end 
suffer  during  that  time  about  10^  collisions.  Direct  experimental  evi¬ 


dence  for  the  existence  of  the  heating  zone  is  furnished  by  observations 

24 

of  Van  de  Poll  and  Westerdijk  \  Figure  5  shows  a  simultaneous  striae  - 
and  normal  photograph  of  a  burner  flame  by  these  authors.  It  can  be 


clearly  seen  that  a  constant  distance  is  maintained  over  the  whole  flame 
surface  between  the  striae  image,  (inner  trace),  which  indicates  tho  be¬ 
ginning  of  the  heating  sone,  and  the  image  of  the  luminous  layer  (outer 
trace),  indicative  of  the  reaction  zone. 

Equations  (30)  end  (31)  account  for  the  observed  dependence  of  the 
flame  front  thickness  on  pressure^.  Since  X,  Cp  and  depend  very 
little  on  pressure,  d{  and  dr  oust  be  approximately  inversely  proportional 
to  P  in  accordance  with  the  experimental  reeulte. 


25 

In  a  recent  paper  ,  i runk-Kemenetaky  drewu  conclusions  from  the 
thermal  theory  on  the  shape  of  a  flame  near  walls  of  a  tube.  It  has  baen 
mentioned  on  page  12  that  according  to  ths  nawar  thermal  theory,  the 
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FIG.  5 

SIMULTANEOUS  NORMAL  AND  STRIAE 
PHOTOGRAPH  OF  A  BUNSEN  FLAME. 
(AFTER  VAN  OE  POLL  AND  WESTERDI JK 

REF.24  ) 


FIG.  8 

FLAME  PROPAGATING  IN  A  CLOSED  CYLINDER 
(AFTER  ELLIS  AND  WHEELER,  REF  26) 


temperature  in  the  reaction  zone  can  be  lower  than  "If  only  by  the  smell 
amount  E  Frank-Komenetsky  concludes  therefore  that  the  flame  front 
con  reach  the  wall  only  with  practically  constant  temperature  Tf  or 
break  up  before  reaching  it,  A  rigorous  theory  would  have  to  he  based  on 
the  two-dimensional  heat  conduction  and  flow  problem  at  the  wall.  Instead, 
as  a  rough  approximation,  Frank-Kamenetsky  concludes  that  no  heat  flow 
parallel  to  the  flame  front  can  take  place  in  the  reaction  zone.  The 
following  cases  are  discussed  in  hia  paper: 

10  The  heat  flow  Cli  ~  Cp  J)  *  Sf  (Tf  —  T0)  perpendicular  to  the 
flame  front  is  greater  than  the  heat  flow  through  the  wall 


A*  =  P  (Tf  ”To)  * 

where  is  the  heat  transfer  coefficient,  end  it  is  assumed  that 
the  outside  of  the  tube  is  kept  at  the  initial  temperature  Te  .  The 
flame  will  then  reach  the  wall  at  such  an  angle  that  the  component  of 
Aj  perpendicular  to  the  wall  equals  A£  ,  Fig  6a  The  er.gle 
must  then  fulfill  the  condition 


sin 


P 


(33) 


^  Cp/0  St 

where  is  taken  positive  if  it  contains  the  unburr. -  gas.  It  is  seer 

that  must  always  be  negative  if  teat  flows  tlirough  the  wall  from 

Inside  to  outside,  Experimenttl  evidence  for  this  behavior  of  a 
flame  is  given,  for  instance  in  a  series  of  flame  photographs  made  by 
Ellis  and  Wheeler  ,  Figure  8;  it  shows  that  the  flsme,  propagating 
in  a  closed  tube,  curves  always  backwards  at  the  wall,  no  matter  how 
the  gas  flow  distorts  Its  shape  during  the  travel. 

-  27  - 


2o  If  heat  flows  from  the  wall  into  the  flame,  e  g.  due  to  catalytic  re¬ 
action  or  by  heating  it  sufficiently  from  outside,  ^  assumes  positive 
values,  Figure  6b, 

3o  The  wall  stays  cold  inside,  (j0»  cp^  Sf)  ;  then  the  flame  cannot 
reach  it  but  will  assume  a  shape  as  in  Figure  7a..  The  distance  d 

will  be  roughly  equal  to  the  thickness  of  the  heating  zone  given 
by  equation  (30),  In  the  case  of  a  flame  burning  outside  a  tube, 
(Bunsen  burner),  the  same  distance  must  appear,  Figure  7b,  The  views 
and  experimental  results  reported  in  a  paper  by  Garside,  Forsyth  and 
Townend^,  particularly  the  analogous  behavior  of  "dead  space"  and 
rounding  of  the  flame  tip  in  burner  flames,  and  of  limiting  diameter 
in  tube  propagation,  are  in  agreement  with  Frank-Kamenetsky '  s  conclu¬ 
sions, 

F o  The  theory  of  Lewis  and  von  Elbe 

It  ha3  been  pointed  out  by  various  authors^ that  the  theory  of 
flame  propagation  must  be  revised,  taking  into  account  the  kinetics  of 
chain  reactions,  and  the  presence  of  active  centers  (free  atoms  and  radi¬ 
cals)  in  the  reaction  zone.  The  first  attempt  in  this  direction  is  due  to 

2,27 

Lewis  and  von  Elbe  ,  in  their  treatment  of  the  decomposition  of  ozone. 
These  authors  recognized  the  importance  of  diffusion  in  flame  proptq aticn, 
and  consequently  based  their  theory  on  the  diffusion  equation 

dN0  j 

Dl^'_d7  (N°*S*)  +  WI°  »*> 
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SHAPES  OF  A  FLAME  NEAR  A  WALL 
(AFTER  FRANK-KAMENETSKY  REF  25) 


for  the  concentration  Ng^  0f  oxygen.  Here  0  is  the  ail  fusion  co¬ 
efficient*  $x  the  flew  velocity  at  X  given  a ppr ox irately  by 


Sx  -  S ♦  — 


T 
1  Ta 


(3t>) 


and  W  the  reaction  rate  of  the  reaction  0  +  03  ~ 1 ^20^,  (a) 


dN 


W  = 


d  t 


°1 _ =  27.e"^ 


(36) 


Since  02  is  present  in  considerable  exco.os,  its  concentration  can  be 
approximately  written 


%  =  N°2  «)  ■  T 


(37) 


The  following  additional  assumptions  are  made; 
(1)  The  reaction 


o3  sr  of  +  o  •  •  -  (M 


is  very  fast  in  comparison  to  (a), 
linked  by  an  equilibrium  equation: 


NoNo 


No, 


*  -  K  -  3.03  .10® 


so  thet  the  concentrations  can  be 


Uml  Jl  e-|2500/T  ,38) 


The  expression  for  the  equilibrium  constant  was  derived  from  data 
by  Kassel. 

(2)  The  sum  of  thermal  and  chemical  energy  is  constant  throughout  the  re¬ 
action  zone.  This  assumption  is  introduced  in  the  form 


where 


e0}  4  N0E0  =  (N0j  4  1.5  N0j  4  0.5 Nq)  cp  (Tf-T) 

is  the  molar  heat  >f  decomposition  c f  ozone  and 


(39) 

the 
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molar  heat  of  recombination  of  oxygen  atoms,  both  at  the  temperature 
T,  and  Cp  the  mean  molar  heat  capacity  of  oxygen  between  T  and  Tf  . 


(3)  Finally,  the  following  equation  is  introduced: 


StN 


°s  1°  J 


(40) 


T. 

which  expresses,  that  all  the  ozone  entering  the  reaction  zone  in 
the  unit  of  time  must  be  decomposed  after  having  passed  through  it. 


For  the  collision  number  2  and  the  diffusion  coefficient  D  expres¬ 
sions  given  by  the  kinetic  theory  of  gases  are  introduced.  The  final  re¬ 
sult  is  given  in  the  form: 


St  =  1.05  X  lo"8Cp(l  +  m)  To  p'/zx^  J 


t. 


- 19300/. 


( 


1.76x10  e 


14  -I2300/t 


+  3.4  x  10  p 


4-7 dTj  (U) 

\  nl  ' 


m  is  the  initial  molar  ratio  of  oxygen  to  ozone.  A  comparison  of 
values  obtained  by  graphical  integration  of  (41)  and  of  experimentally 
determined  transformation  velocities  is  given  in  Table  1. 


TABLE  1  (after  Lewis  and  von  Elbe) 


p 

To 

Tf 

m 

cole. 

txper. 

mm 

#K 

#K 

cm^sec 

cm/t9c 

624 

300 

1239 

3.054 

253 

55 

2560 

427 

1343 

3-054 

451 

158 

595 

302 

1922 

1.016 

333 

160 

3760 

468 

2044 

1.016 

664 

747 
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FIG.  9 

STRUCTURE  OF  REACTION  ZONE  OF  OZONE  FLAME 
49.6%  OZONE  IN  OXYGEN,  P  =  3760MM,  T0»468°K, 

Tf  =  2044°K  (AFTER  LEWIS  AND  V.  ELBE, REF  2  AND  27) 


For  the  last  experiment  listed  in  the  table,  Lewi.1  and  von  Elbe  give 
a  graphical  representation  of  the  structure  of  the  reaction  zone,  Figure  9. 
It  seems  rather  improbable,  however,  that  the  reaction  rate  should  reach 
its  maximum  at  so  low  a  value  of  the  temperature,  as  shown  in  this  graph. 

The  theory  of  Lewis  and  von  Elbe  constitutes  the  first  successful  cal¬ 
culation  of  the  absolute  value  of  the  transformation  velocity  S^.  The 
fair  agreement  between  experimental  and  theoretical  results  could  be 
achieved  because  the  particular  kind  of  reaction  selected  by  the  authors 
allowed  the  introduction  of  simplifying  assumptions.  The  fundamental 
assumption  of  this  theory,  namely,  the  constancy  of  the  sum  of  thermal  and 
chemical  energies,  was  later  shown  by  Zeldovich  and  Frank-Kamenetsky^  to 
hold  under  the  condition  that  the  coefficients  of  temperature  conduction 
and  diffusion  are  equal  (see  p.  16  )e  They  showed  further  the  equivalence 
of  the  equations  of  diffusion  and  of  heat  conduction  under  the  same  condi¬ 
tion.  Both  the  theory  of  Lewis  and  von  Elbe  and  the  thermal  theory  are 
based  in  addition  on  a  kinetic  law  of  the  formW*S  e  Therefore,  the 

two  theories  should  give  similar  results  and  the  theory  of  Lewis  and  von 
Elbe  constitutes  a  particular  case  of  the  more  general  thermal  theory  of 
Zeldovich. 
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Go  Formulation  of  the  thermal  theory  for  more  complex  reaction  mechanism 
In  a  qualitative  way,  Zeldovich  and  Semenov'  have  discussed  the  prob¬ 
lem  of  dealing  with  chain  reactions  in  flame  propagation  They  distinguish 
two  main  types  of  reactions! 

Type  li  The  life  tine  of  the  active  centers  is  short  compared  with  the 

duration  of  the  reaction  In  this  case,  the  instantaneous  concen¬ 
trations  of  the  active  centers  can  be  determined  with  sufficient 
accuracy  by  equilibrium  with  the  reacting  substances  and  are 
tiniquely  connected  with  the  concentrations  of  these  substances. 

The  reaction  rate  can  be  expressed  in  this  case  by  a  "classical" 

“TfT 

law  of  the  formW  =  Sg  ,  in  terms  of  the  concentrations  of  the 
reacting  substances  alone.  For  this  type  of  reaction,  the 
theories  of  Lewis  and  von  Elbe  and  of  Zeldovich  and  Frank-Kamenetsky 
can  be  applied., 

Type  2:  The  life  time  of  the  active  centers  is  of  the  same  order  as  the 
time  of  reaction  or  larger.,  This  condition  leads  to  the  auto- 
catalytic  type  of  chain  reaction,  This  class,  with  the  combustion 
of  tydrogen  and  the  frydrocarbone  as  outstanding  examples,  presents 
difficulties  for  the  theory  of  flame  propagation  which  have  not  yet 
bren  overcome.  The  complications  are  due  only  in  part  to  the 
purely  mathematical  side  of  the  problem,  to  a  greater  part,  however, 
to  the  lack  of  definite  knowledge  about  the  mechanisms  of  the 
reactions. 
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tormulsted  quite  generally,  as  has  been  done  for  instance  by  Damkohlei  , 
(cf.  ref,  9  pp-  115-117) r  the  problem  consists  in  finding  a  simuitenecus 
solution  of  the  differential  equation  of  heat  conduction,  equation  (6), 
and  of  as  many  diffusion  equations  as  there  are  different  kinds  of  par 
tides  taking  part  in  the  reactions*  In  these  equations  appear  the  terms 
which  give  the  generation  of  teat  and  the  creation  or  destruction  of  par¬ 
ticles,  which  depend  expliciteiy  on  the  temperature  and  the  concentrations, 

not  on  the  independent  variable  X  .  This  has  led  Zeldovich  and  Jrnnk- 
19 

Kamenetsky  to  eliminate  X  entirely  from  the  equations,  by  introducing 
the  following  set  of  variables?  as  independent  variables,  the  non-dimensional 
temperature  8  and  the  non-dimensional  concentration  V; (for  instance  by 

T  -Ta  n  -  no 

netting  0  =  — — —  ,  V  - - a- 

VTo  nf“no 


j _ a _ a. 


X-* 


de 

dx 


*  s  dx 


(42) 


where  is  a  constant  of  the  dimension  of  a  length  (e<.g0  equal  to  the 

jr  ^ 

thickness  of  the  heating  zone  T  =  —  = - ).  In  these  variables,  the 

st  cp^st 

differential  equations  of  heat  conduction  and  diffusion  become 


Tj  ^  -  T|^  +  X^O.U)- JX'-J’t©)  =0  (43) 

K  X  +  xsf(e,u)  --o  (t4> 


*  These  equations  lack  generality,  however,  in  that  the  factor  of  the  second 
term  cannot  be  unity  in  all  equations  if  the  condition  D  =  X  is  not  ful¬ 
filled,  and  in  that  the  same  function  1(6, V)  need  not  apply  for  the 
generation  of  beat  and  of  particles. 


- 
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where  Xj.Xg  and  JA  are  non-dimensional  parameters,  containing  emcng 
other  constants  the  transformation  velocity  S^,fis  a  dimensionless 
function  expressing  the  dependence  of  the  reaction  rate  on  temperature 
and  concentration,  and  ^  (0)  expressing  the  heat  loss  due  to  walls  and 
radiation.  In  the  general  case,  several  equations  (44)  have  to  be  set 
up.  With  the  indicated  choice  of  0  and  V  the  simple  boundary  condi¬ 
tions 


11=0,  t  =  °  <«'  e  =v  =  o  and  e  =  v  =  i  (45) 

are  obtained.  The  authors  give  explicit  solutions  only  for  the  following 
cases: 

For  reactions  characterised  above  under  type  1,  ond  neglecting  her*  'loss, 

4  n  1  q 

there  results  the  thermal  theory  *  *  »  leading  to  the  general  solution 

St = \/2p,f^h  (26) 

which  has  been  discussed  on  page  18 „ 

For  the  other  extreme,  an  isothermal  autocatalytic  reaction  ("cool 
flame", diffusional  propagation)  only  equation  (44)  in  the  form 

has  to  be  solved.  The  authors  give  the  solutions 
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St=  yf /f3  =2/f*  (46) 

for  a  kinetic  law  w  =(^>  "1/  (l~V)  (46a) 

where  v^>  is  the  coefficient  of  auto-accelereticn,  and 

St  Ml -20 (47) 

for  a  kinetic  law 


w  =  ^  UZ  (1  — V)  “  a^)V(l~V)  (47e) 

corresponding  to  auto-acceleration  cf  second  order  with  destruction  of 
active  centers.  a  is  the  ratio  of  rate  of  destruction  to  rate  of  re¬ 
action.  It  is  seen  that  in  the  latter  case  propagation  is  limited  to  the 
range  a  <  -y  .  it  is  pointed  out  by  the  authors,  tliat  kinetic  conditions 
in  ell  real  cases  are  too  c-omplicuteo  to  allow  an  analytical  solution  of 
equations  (43),  (44) 


H.  Hydrodynamic  aspects  of  flame  propagation 

The  theory  as  discussed  up  to  here  neglectE  the  possible  effects  of 
motion  in  the  gas  Gas  flow,  however,  is  always  present,  either  due  only 
to  the  combustion  process  itself,  or  in  addition  introduced  by  external 
conditions. 

The  Hydrodynamics 1  treatment  of  combustion  under  the  assumption  of 
strictly  one-dimensicnal  flow  leads  to  the  Hugoniot  curve,  Fig.  10,  which 
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gives  the  locus  of  the  final  states  for  a  given  initial  state,  in  terms  of 

pressure  p  and  specific  volume  V.  A  detailed  review  of  this  theory  has 

43 

been  given  in  a  recent  report  so  that  the  discussion  here  can  be  brief. 

Slow  combustion  (deflagration)  corresponds  to  the  portion  AB  of  the 
Hugoniot  curve.  The  velocities  of  the  unburnt  and  Vg  of  the  burnt 
gas  relative  to  the  flame  front  are  both  zero  in  point  A.  In  going  from 
A  to  B,  both  increase;  in  point  B  Vt  reaches  a  maximum  but  remains 
always  smaller  than  the  velocity  of  sound  in  the  unbumt  gas,  Cj  and  Vg 
becomes  equal  to  the  velocity  of  sound  in  the  burnt  gas,  C2  .  It  can  be 
shown  that  the  flow  states  corresponding  to  points  beyond  B  ( Vl<CJ  ,V2>cg) 
cannot  be  realized.  Which  point  between  A  and  B  will  correspond  to  the 
actual  combustion  process  is  not  determined  by  hydrodynamics.  In  the  case 
of  strictly  one-dimensional  flow  with  a  plane  flame  front  Vt  is  identical 
with  the  transformation  velocity  .  From  the  point  of  view  of  the 
thermal  theory  of  flame  propagation  it  has  been  always  assumed  that  St  *  V| 
and  Vg  are  so  small  that  the  combustion  process  can  be  represented  by  a 
point  very  near  to  A  on  the  Hugoniot  curve.  In  other  words  it  has  been 
assumed  that  Vx  and  Vg  C2  .  This  assumption  is  not  too  re¬ 

strictive  in  many  cases,  since  the  observed  values  of  are  very  small 
compared  with  the  velocity  of  sound  in  the  unburnt  gas.  As  long  as 
VT  ^  ^  Cj  and  Vg  ^  ^  Cg  ,  the  pressure  gradient  across  the  flame  front 
will  be  very  small  (ref.  2,  p.  163)  and  its  effect  on  heat  conduction, 
diffusion  and  chemical  kinetics  can  be  neglected. 
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For  larger  values  of  Vj  and  Vg  an  influence  of  the  flow  on  the 
transformation  velocity  cannot  be  excluded  any  more.  No  theory  on  this 
effect  exists  yet.  It  appears  reasonable  to  assume,  however,  that  the 
rise  of  pressure  and  temperature  due  to  the  compression  wave  traveling 
ahead  of  the  flame  front  will  lead  to  an  increased  transformation  velocity P 
This  acceleration  can  take  place  only  until  point  B  is  reachedc  Propaga¬ 
tion  processes  with  still  higher  velocities  correspond  to  the  portion  CD 
and  beyond  D  of  the  Hugoniot  curve  (detonation).  It  has  been  shoim^6 
that  under  the  assumption  that  the  detonation  process  consists  of  a  shock 
wave  followed  by  combustion,  and  assuming  strictly  one-dimensional  flow, 


stable  propagation  corresponds  only  to  point  D  ,  where  Vj  is  a  minimum, 

and  Vj  ^  Cj  ,  Vg=  c2  »  This  state  is  reached  only  through  unstable 

transition  states.  Detonation  will  not  be  discussed  further  in  this  re- 

2  9,9a 

port  and  reference  is  made  to  the  texts  of  Lewis  and  von  Elbe  and  Jost 
for  reviews  of  this  subject-. 

The  assumption  of  one-dimensional  flow  is  actually  never  fulfilled 
for  a  flame  propagating  in  a  pipe.  Jost  (ref.  9,  p.  92-96)  has  stressed 
the  fact  t'aat  the  unburnt  gas  must  be  in  a  vortex-like  motion  in  order  to 
enable  the  front  of  a  flame  propagating  in  a  tube  to  assume  the  convex 
shape  which  is  always  observed.  (This  shape  will  be  due,  at  least  in  part, 
to  the  condition  imposed  on  the  flame  near  the  wall,  p.  26  ).  A  hydrody- 
namical  treatment  of  this  complicated  motion  has  not  yet  been  attempted. 

The  effect  of  this  motion  on  the  flame  speed  is  accounted  for  theoretically 
only  In  an  increase  of  flame  surface  without  a  change  in  transformation 
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velocity c  However,  the  motion  of  the  gaB  changes  the  shape  of  the  flame 
front  so  profoundly  that  it  becomes  difficult  or  even  impossible  to  cal¬ 
culate  the  transformation  velocity  from  the  observed  propagation  speed » 
This  has  precluded  a  check  of  the  validity  of  the  theory  in  this  case. 


A  large  effect  of  flow  on  flame  propagation  is  finally  to  be  expected 

where  turbulent  motion  is  present.  That  such  turbulence  may  even  originate 

in  the  flame  itself,  due  to  a  peculiar  instability,  has  been  demonstrated 
29 

by  Landau  ,  His  treatment  is  based  on  the  following  assumptions: 

(a)  The  dimensions  which  characterize  the  tydroaynamicel  problem  (such  as 
tube  diameter,  dimensions  of  bodies  in  contact  with  the  gas)  are  large 
compared  with  the  thickness  of  the  combustion  layer,  bo  that  the  latter 
can  be  treated  as  a  surface  of  discontinuity, 

(b)  The  propagation  velocity  is  small  compared  with  the  velocity  of  sound, 
allowing  the  gas  to  be  treated  as  incompressible, 

(c)  Viscosity  is  neglected,  restricting  the  validity  of  the  result  to 

a  v,  ov, 

large  values  of  the  Reynolds  numbers  — —  ,  — —  where  0  is  a  char- 

U  i  V  z 

acteristic  dimension,  V,  and  V2  sure  the  velocities  of  the  unburnt 
and  burnt  gas,  and  Vj  and  '\J  are  the  kinematic  viscosities  of  un¬ 
burnt  and  burnt  gas. 


As  is  usual  in  the  theory  of  flame  propagation,  the  origin  of  the  co¬ 
ordinate  system  is  fixed  in  a  small  portion  of  the  combustion  layer,  which 
is  considered  as  a  plane  discontinuity  surface;  the  X  -axis  points  in 
the  direction  of  the  unperturbed  gas  movement,  which  is  assumed  to  be 
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velocity*  However,  the  motion  of  the  gas  changes  the  shape  of  the  flame 
front  so  profoundly  that  it  becomes  difficult  or  even  impossible  to  cal¬ 
culate  the  transformation  velocity  from  the  observed  propagation  speed. 
This  has  precluded  a  check  of  the  validity  of  the  theory  in  this  case. 


A  large  effect  of  flow  on  flame  propagation  is  finally  to  be  expected 
where  turbulent  motion  is  present.  That  such  turbulence  may  even  originate 
in  the  flame  itself,  due  to  a  peculiar  instability,  has  been  demonstrated 
by  Landau  ,  His  treatment  is  based  on  the  following  assumptions: 

(a)  The  dimensions  which  characterize  the  hydrodynaroicel  problem  (such  as 
tube  diameter,  dimensions  of  bodies  in  contact  with  the  gas)  are  large 
compared  with  the  thickness  of  the  combustion  layer,  so  that  the  latter 
can  be  treated  as  a  surface  of  discontinuity, 

(b)  The  propagation  velocity  is  small  compared  with  the  velocity  of  sound, 


allowing  the  gas  to  be  treated  as  incompressible, 

(c)  Viscosity  is  neglected,  restricting  the  validity  of  the  result  to 

a  V,  a  V2 

large  values  of  the  Reynolds  numbers  -r: —  ,  -rr—  where  a  is  a  char- 

l/i  V  2 

acteristic  dimension,  V,  and  Vg  are  the  velocities  of  the  unburnt 


and  burnt  gas,  and  Vj  and  '[)  are  the  kinematic  viscosities  of  un¬ 
burnt  and  burnt  gas. 


As  is  usual  in  the  theory  of  flame  propagation,  the  origin  of  the  co¬ 
ordinate  system  is  fixed  in  a  small  portion  of  the  combustion  layer,  which 
is  considered  as  a  plane  discontinuity  surface;  the  X  -axis  points  in 
the  direction  of  the  unperturbed  pas  movement,  which  is  assumed  to  be 
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perpendicular  to  the  combustion  layer,,  The  index  1  refers  to  the  unfcurnt 
gas(X^O),  the  index  2  to  the  burnt  gas(X^0K  A  snail  perturbation 

i 

velocity  V  is  superposed  over  the  unperturbed  movement  V  .  The  com- 
/ 

ponents  of  V  must  fulfill  the  differential  equations  (neglecting  small 
terms  of  second  order) 


(continul ty) 


(48) 


Qvy  1  3p 

- *.  4-  v  - Cl  z - - - 

3x  f  3x  (Euler)  (49) 

3y;  i  V 

3  t  3  X  ^r,  5  Y 

’.v.ere  J>  is  the  density  and  p  the  perturbation  pressure-  Equations 
(48)  and  (49)  lead  to 


(50) 


Solutions  are  sought  which  are  periodic  in  time  and  in  Y  ,  by  setting 

Iky  +nt 

them  proportional  to  C  „  Instability  will  then  correspond  to 

the  case  R.  (n»  o  .  The  solutions  fulfilling  equations  (48) ,  (49) 
and  (50)  are 
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J 


V,, *  A« 


iky  +  kx  +■  At 


Vly  siAe 


Pi 


iky  ♦  kx  +  At 


' s  A/i  Hr+v')e 


iky  +  kx  +  A  t 


(51) 


In  the  unbumt  gas,  and 


vzY. 


*2  Y 


0  e iky  -  k*  +  At  +  c#  iky  +  At  -  X 


iky  -  Kx  4- At  r  iA  iky  ♦  At 

iBe  - C  — -  e 

kvt 

A  \  iky  -  kt  +  At 


A 


(51a) 


K  -  -  BA  (v*  -  — ) * 


for  the  burnt  gaeo  It  will  be  seen  that  the  solutions  are  set  up  so 
that  they  vanish  for  jxj  - ►  oo  ,  if  R*  (k)  0  ,  and  R»  (A)^O. 

This  is  of  course  necessary  if  they  should  indicate  that  a  distur¬ 
bance  small  everywhere  else  becomes  large  within  the  combustion  zone . 

For  this  reason  the  particular  solution  with  the  factor  C  ,  obtained  by 
setting  p*  s  0  cannot  be  used  for  the  unburnt  gas. 

At  the  discontinuity  (X=0)  the  following  boundary  conditions  must 
be  fulfilled:  based  on  assumptiou  (b)  it  is  assumed  that  the  propagation 
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velocity  is  not  affected  by  the  disturbance,  The  displacement  of  the 
discontinuity  in  the  X  -direction 

.  .  .  _  iky  ♦  At 

^  ( Y,t )  =  De 

is  therefore  connected  with  the  disturbance  velocities  by 


viX  s  V£X 


for  X  *0 


(53) 


To  the  same  approximation)  the  pressure  is  assumed  to  be  continuous 
across  the  discontinuity: 


P  Z  Q  for  X  =  0 
1  2 


(54) 


Finally,  the  tangential  component  of  the  velocity  must  be  continuous  at 
the  discontinuity o  The  angle  of  inclination  of  the  discontinuity, 

Ot 

assumed  to  be  small  in  the  initial  stage  of  perturbation, is  ■■  Y  i  , 

OY 

Thus,  to  a  first  approximation,  the  tangential  velocity  is  composed  of  the 


|  * 

projection  of  V  equal  to  Vy  ,  and  the  projection  of  V  equal  to  V 


This  gives 


vlY  +  V«  = 

dY 


vt  Y  *  V2  “o  f0r  X  S  0 

dv 


as  last  boundary  condition. 


(55) 


The  conditions  (53) »  (54)  anu  (55)  lead  to  fou\*  linear  homogeneous 
equations  for  the  constants  A,  B,  C,  D,  w!  .ch  tev*>  a  non-trivial 
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solution  only  if  the  equation  for  A  : 


il*  ( v,  +  v8 )  +  2  A  v,  v£  k  +  k* v,  vt  (v,  -  Vt)  sO  (56) 

is  fulfilled,  (where  the  continuity  condition,  J>{  V(  *J>  V  ,  has  been 
taken  into  account).  In  the  case  of  combustion,  there  is  always 
^  i  Vj  ^  V£  .  Equation  (56)  has  then  always  solutions 
for  which  Re  (A)  ^  0  (considering  Rs  (k)  ^  0).  Unstability  under 
the  assumed  conditions  is  thus  proved.  Landau's  conclusions  from  this 
result  are  as  follows: 


The  "thermal"  regime  of  combustion,  on  which  the  usual  theories  of 

* 

flame  propagation  are  based  ,  is  valid  only  inder  conditions  which  de¬ 
pend  essentially  on  viscosity  (small  Reynolds  number,  contradictory  to 

the  assumptions  made  in  the  derivations  of  the  instability).  Under 

** 

conditions  of  large  Reynolds  numbers,  however  ,  instability  will  lead 
to  a  diffuse  combustion  zone  in  which  turbulence  will  cause  heat  trans¬ 
fer  considerably  in  excess  of  that  due  to  conduction.  Since  turbulent 
motion  depends  essentially  on  the  characteristic  dimensions  of  the 

*  this  obviously  includes  "diffusional"  propagation  p.  34 

*»  Though  not  stated  clearly  in  Landau's  paper,  his  implication  is  ob¬ 
viously  that  turbulence  need  not  exist  in  the  gas  before  the  flame 
has  passed.  The  unbumt  gas  may  also  be  at  rest  relative  to  the 
tube.  Turbulence  will  appear  always  in  the  flame  zone  if  the 
Reynolds  numbers  ^  t  are  large,  where  v,  ,  v2  are  the 

V,  ’  V2 

gas  velocities  relative  to  the  combustion  zone. 


problem,  Landau  concludes  that  in  this  "convective"  regime  of  combustion 
the  width  of  the  combustion  zone  and  the  propagation  velocity  will  de¬ 
pend  on  the  characteristic  dimensions.  More  explicitly,  for  flame  propaga¬ 
tion  in  a  tube,  and  under  conditions  where  viscosity  and  heat  conduction 
(or  diffusion)  cannot  have  influence,  Landau  concludes  from  dimensional  con¬ 
siderations  that  the  thickness  d  of  the  flame  zone  must  be  proportional  to 
the  tube  diameter  (5  and  the  propagation  velocity  must  be  proportional  10 
where  T  is  a  reaction  time  dependent  only  on  the  chemical  kinetics. 


This  result: 


d 


v  ~  4- 

t 


(57) 


Is  said  to  be  valid  for  V  small  compared  to  the  velocity  of  sound  and 
for  large  Reynolds  numbers. 


Though  this  result  will  have  to  be  accepted  with  some  caution,  exper¬ 
imental  data  confirm  the  increase  of  propagation  velocity  with  tube  di- 

30 

ameter.  Figure  11  shows  the  results  of  measurements  by  Coward  and  Hartwell 
(ref.  9 ,  p.  99).  Two  distinct  regions  can  be  observed  in  the  graph:  the 
increase  of  observed  flame  speed  up  to  a  diameter  of  about  10  cm  could  be 
explained  by  the  authors  by  determining  the  surface  of  the  flame  front  from 
flame  photographs,  Figure  12.  Table  2  shows  that  tbe  true  transformation 
velocity,  obtained  by  multiplying  the  observed  speed  with  the  ratio  of  tube 
cross  section  to  flame  surface,  is  constant. 
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TABLE  2 


Flame  speeds  in  10$  Methane-air  mixture  in  tubes  of  various  diameters 

(Coward  and  Hartwell,  ref.  30) 


TUBE 

DIAMETER 

FLAME 

MOVEMENT 

OBSERVED  FLAME 
SPEED 

FLAME 

SURFACE 

TRANSFORMATION 

VELOCITY 

cm 

ciVsec, 

2 

cm 

cm/sec , 

10 

horizontal 

111 

300 

29 

10 

n 

71 

189 

29 

5 

n 

92 

66 

27 

5 

n 

61,5 

48  5 

25 

2o5 

n 

?1,5 

12  .6 

28 

2o5 

n 

63 

11,0 

28 

2,5 

M 

59 

10,4 

28 

5 

Upward 

68 

48 

28 

5 

ii 

92,5 

66,5 

27 

5 

Downward 

61 

46 

26 

5 

n 

58 

63  5 

26 

This  range  corresponds  to  the  regime  called  "thermal"  by  Landau,  For 
diameters  between  10  and  20  cm,  the  graph  Figure  11  indicates  an  unstable 
transition  and  for  still  larger  diameters  an  approximately  linear  increase 
of  propagation  velocity.  Coward  and  Hartwell  point  out  that  in  this  range 
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FIG.  II 

FLAME  SPEEOS  IN  METHANE-AIR  MIXTURES  IN 
TUBES  OF  VARIOUS  OIAMETERS  (AFTER  COWARD 
AND  HARTWELL  REF  30) 
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FIG  12 

PHOTOGRAPHS  OF  METHANE-AIR  FLAMES  PROPAGATING 
IN  TUBES  OF  VARIOUS  DIAMETERS. 
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it  was  impossible  to  estimate  the  flame  surface,  since  it  became  irreg¬ 
ular  and  diffuse.  The  beginning  of  this  irregular  appearance  can  al¬ 
ready  be  observed  in  the  photograph  for  diameter  10  cm,  Figure  12.  It 
seems  reasonable  to  assume  that  this  range  is  identical  with  Landau's 
"convective"  regime. 

Landau  includes  in  his  paper  a  discussion  of  the  limitations  im¬ 
posed  on  the  propagation  velocity  by  the  laws  of  gasdynamic3  (cf.  p. 35) 

His  proof  of  the  instability  of  processes  corresponding  to  points  beyond 
Bon  the  Hugoniot  curve  (Figure  10)  appears  to  be  unnecessarily  compli¬ 
cated.  The  subject  is  treated  in  a  more  satisfactory  way  in  ref.  43  and 
43a.  It  should  be  borne  in  mind  that  the  limiting  velocity  corresponding 
to  point  B  lies  already  outside  the  range  of  validity  of  the  theory  of 
Landau,  ns  the  velocity  of  the  burnt  gas  Vg  relative  to  the  flame  front 
becomes  equal  to  the  velocity  of  round  at  this  limit.  Therefore,  de¬ 
viations  from  the  simple  expressions  (57)  should  be  expected  already  for 
lower  velocities.  In  u  qualitative  way  Landau' s  theory  explains  many  of 
the  difficulties  experienced  by  workers  in  the  field  of  flame  propagation 
in  tubes  (ref.  9,  page  101),  and  the  instability  postulated  by  Landau 

mey  be  responsible  1'or  the  increuse  in  flame  velocity  obtained  by  intro- 

31  32 

ducing  constrictions  or  wire  spirals  **  into  the  tubes. 

The  influence  of  turbulence  on  the  velocity  of  flame  propagation  has 

„  ,  13,13a  ,  33,33a 

been  investigated  by  Damkohler  and  Lhchelxin  .  Both  authors 

base  their  treatment  on  the  description  of  turbulent  flow  by  means  of  a 
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characteristic  length  41  an  average  fluctuation  velocity  V  ,  and  an 

n 

exchange  quantity  €  <>  Damkohler  defines  these  quantities  according  to 

the  phenomenological  theory  of  Prandtl  (cf.  ref.  34):  the  mixing  length^ 
is  defined  by 


(53) 


where  I  V  I  is  the  time  average  of  the  absolute  value  of  the  instantaneous 
fluctuation  velocity  V,  and  the  absolute  value  of  the  gradient  of 

dY 

the  flow  velocity  U  in  a  direction  Y  perpendicular  to  the  flow.  The 
exchange  parameter  6  was  introduced  by  Prandtl  in  analogy  to  the 
kinetic  theory  of  gases:  in  a  quiescent  gas  kinematic  viscosity  U  ,  co¬ 
efficient  of  diffusion  D  and  coefficient  of  temperature  conduction  X 
are,  to  a  first  approximation,  equal  end  given  by  (cf.  p.  16  and  ref.  22) 


U  =  0  =  X=i  He  (59) 

whare  Si  is  the  mean  free  path  and  C  the  mean  molecular  velocity. 
Prandtl  showed  that  in  a  turbulent  medium  these  coefficients  heve  to  be 
replaced  by  the  quantity  6  ,  which  is  defined  in  an  analogous  way  by  6 

6  *  lY'l  (60) 

Shchelkin  adopts  definitions  in  accordance  with  the  statistical 
theory  of  Taylor  and  Karman  (cf.  ref.  35):  as  characteristic  length 
he  introduces  the  scale,  defined  by 


where  RW  is  the  correlation  coefficient 


R  * 


(61a) 


between  the  fluctuation  velocities  V*  and  vj  at  two  points  a  distance 
X  apart.  The  average  fluctuation  velocity  is  defined  by  Shchelkin  as 
the  R.M.S.  value  of  V  .  Unfortunately  Shchelkin  makes  no  distinction 
between  scale  L  and  mixing  length  JL  in  his  treatment,  using  expression 
(60)  in  spite  of  the  fact  that  in  the  statistical  theory  the  coefficient 
of  exchange  €  is  not  related  to  the  scale  by  this  formula.  Because  both 
Daukohler's  and  Shchelkin' s  theories  are  of  a  rather  crude  and  qualitative 
nature,  these  discrepancies  of  definitions  will  be  disregarded  in  the 
following  discussion. 


n 

The  following  cases  are  treated  by  Damkohler  and  by  Shchelkin: 

I 

(a)  The  mixing  length  X  is  small  compared  with  the  thickness  of  the 
flame  front  d  ,  JL«d 


In  this  case,  the  shape  of  the  flame  front  will  not  be  affected  by 
turbulent  motion.  The  only  effect  will  be  an  increase  of  all  transport 
processes,  like  heat  conduction  and  diffusion.  It  follows  from  dimen¬ 
sional  analysis  that  the  transformation  velocity  is  proportional  to 
\  /  X  or  /  D  ,  where  7“'  is  a  characteristic  time  . 

V  f  V  f 


*  The  expressions  for  Sj  of  the  older  thermal  theory,  pages  7&8,  as  well 
as  those  of  the  new  thermal  theory,  formula  (26)  and  those  for  "diffu¬ 
sions!  "  propagation,  (46)  and  (47),  have  all  this  general  form. 
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If  this  time  is  the  reaction  time,  depending  only  on  the 

physicochemical  properties  of  the  gas  mixture.  Since  X  or  D  has  to  be 

* 

replaced  by  6  in  a  turbulent  flow,  Damkohler  writes  for  the  propaga¬ 
tion  velocity  under  turbulent  conditions,  if  JL 


'turb. 


=  S 


lam 


(62) 


where  S(om  is  identical  with  the  transformation  velocity  in  a 
laminar  flow. 


In  order  to  account  for  the  fact  that  the  propagation  velocity  must 
become  equal  to  Sjom  for  zero  turbulence,  Shchelkin  assumes  that  the 
total  coefficient  of  heat  exchange  is  the  sum  of  temperature  conductivity 
X  and  turbulent  exchange  quantity  6  .  This  leads  to  the  alternative 

expression 


Sturb  =  Slom  \fx~  S  Slom'/1  +  "x  (63) 

Shchelkin  points  out,  however,  that  the  case  X  <^<(d  is  not  likely  to 
occur  in  practice,  as  the  width  of  the  flame  front  is  normally  of  the 
order  0.1  mm  and  mixing  lengths  smaller  than  this  are  not  encountered  under 
the  conditions  of  flow  in  practical  applications.  Only  for  slow-burning 
mixtures  the  thickness  of  the  flame  front  can  be  considerably  larger  and 
in  this  case  the  condition  it  could  be  fulfilled  more  easily. 
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I  . 

(b)  The  intermediate  range.  Jt  \  d.  ig  discussed  only  by  Shchelkin. 

In  this  case,  heat  conduction  (and  diffusion)  will  still  be 
governed  by  the  coefficient  of  turbulent  exchange  €  sJL-V  ,  but  the 
time  T  in  the  general  expression  for  the  transformation  velocity 

fx 

v-^p“  will  no  longer  be  equal  to  the  reaction  time  under  laminar 
flow  conditions  Shchelkin  asserts  this  to  be  so  because  in  a  certain 
plane  parallel  to  the  flame  front  and  within  the  combustion  zone  turbu¬ 
lent  mixing  will  bring  together  gas  of  the  whole  temperature  range 

a'  dT  dT 

to  T  +.JL,  „  If  the  temperature  gradient  is  replaced  approxi¬ 
mately  by  ,  it  is  seen  that  the  effect  of  turbulence  on  the 

reaction  rate  will  become  noticeable  if  cannot  be  neglected 

compared  with  T  .  Shchelkin  points  out  that  the  time  of  reaction , 
proportional  to  6  under^laminar  conditions,  will  increase  by  be¬ 
coming  proportional  to  e  R[T^,T'-T*,J  .  (This  expression  cannot 

be  considered  correct  as  it  corresponds  to  the  maximum  reaction  time,  not 
to  the  average.  However,  considering  the  fact  that  the  regions  at  higher 

temperature  T  +  A  (T,  -To)  correspond  to  regions  occupied  by  burnt  gas, 

d 

turbulence  will  cause  some  increase  in  average  reaction  time.)  The 

author  asserts  that  by  increasing  Si  more  and  more,  a  point  will  be 

reached  where  the  time  of  chemical  reaction  will  become  so  large  that  the 

rate-determining  factor  will  be  the  time  of  turbulent  mixing,  of  the  order 
0' 

-=■  instead. 


In  this  way,  Shchelkin  believes  that  a  continuous  change  from  case 
M  ft'«d,  to  case  (c  >,f>>d  ,  can  be  established. 
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(c)  Mixing  length  large  compared  Tilth  (laminar)  thickness  of 
flame  front  d t  &'>><)  . 

Both  authors  assert  that  in  this  case  the  effect  of  turbulent  mo¬ 
tion  will  consist  essentially  only  in  a  distortion  of  the  flame  front 
and  therefore  in  an  increase  of  its  surface  and  of  the  effective  propa¬ 
gation  velocity  in  the  same  proportion.  (Semenov^*  considers  this  to  be 
the  only  possible  effect  of  turbulence  on  flame  propagation.) 

For  small  values  of  V*  Shchelkin  adopts  the  same  simplified  descrip- 

N 

tion  of  the  flame  front  already  used  by  Damkohler,  Figure  13 .  The  com- 

ponent  of  V  in  the  direction  of  flame  travel  is  assumed  to  create  rapidly 

changing  cone-shaped  elevations  and  depressions  in  the  flame  front.  The  area 

a'2 

of  the  base  of  these  cones  is  of  the  order  Ji  ,  and  the  height  of  the 

_ a _ ^  _ x _  _  i  _  .» _ _ _ .  i _  _ 1  i* _  _ i  _ 


order  V 


e —  ,  since  -x-  is  the  average  time  a  uniform  velocity 

^lom  ^lom 


V  is  operative  to  build  up  a  cone.  The  ratio  of  the  cone  surface  Ac  to 
the  base  A^  is  thus 


—  \  2 


_V_ 

S|om, 


where  B  is  a  constant.  Assuming  the  true  propagation  velocity  Sjom  to 
remain  unchanged  by  turbulence,  the  pparent  propagation  velocity 
will  be  increased  by  the  ratio  of  the  surfaces,  that  is 


Sturt 

S|om 


1  tB 


V'  \f 


-  50  - 


For  large  values  of  V  this  typerfcolic  law  goes  over  asymptotics lly 


into  e  linear  relstioi 


^turb^  V 


(b4a) 


which  had  been  alreudy  postulated  by  Lamkohler.  (An  objection  can  be 

f 

raised  against  this  theory:  for  sufficiently  large  values  of  V  ,  the 
flame  front  will  be  disrupted  and  propagation  will  no  longer  be  possible, 

f 

This  upper  limit  may  depend  not  only  on  V  ,  but  also  on  the  mixing 
length  Ji,  „ )  Shchelkin  discusses  the  case  v»  S|Qm  separately.  He 
describes  the  combustion  zone  in  this  cr.se  as  consisting  of  a  number  of 
islands  of  unburnt  mixture  surrounded  by  burnt  gas  and  decreasing  in  size 
while  passing  through  the  combustion  zone,  Figure  14,  (As  indicated  in 
this  figure, there  must  be  a  continuous  flame  front  toward  the  unburnt  gas, 
otherwise  combustion  would  cease,  as  mentioned  above*)  Shchelkin  asserts 
that  under  these  conditions  heat  exchange  is  determined  entirely  by  the 

coefficient  of  turbulent  exchange  6  V  and  the  rate  of  burning  is  de- 

1 

termined  only  by  the  time  of  turbulent,  mixing,  of  the  order  ~  .  In- 

V*  /f 

troduced  into  the  general  expression  for  the  propagation  velocity  , 

this  gives 


^turb  ^ 


k  i  :  7 


that  is  the  same  linear  law  (b4a)  previously  derived.  This  rather  crude 


derivation  can  hardly  be  considered  correct;  a  computation  of  the  com- 

bustion  time  of  a  gas  volume  introduced  by  Shchelkin  to  justify 

^  .  O' 

further  the  identification  of  /  with 


,  is  evidently  incorrect. 


•  I 


1 


I 


From  expression  (64)  Shchelkin  draws  the  following  conclusions: 
lo  If  the  mixing  length  SL  is  large  compared  with  the  thickness  of  the 
flame  front,  the  turbulent  propagation  velocity  Sturjj  depends  only 
on  V  ,  not  on  Jl . 


2.  For  values  of  7  of  the  order  of  S|om  t  the  relative  increase  of 


propagation  velocity 


'turb 


Mom 


is  larger  for  slow-burning 


mixtures,  for  instance  for  mixtures  near  the  limits  of  propagation. 


For  V  large  compared  with  S|om  ,  depends  on  V  accord¬ 

ing  to  the  linear  relation  (64a);  the  factor  of  proportionality^) 
depends  probably  on  the  "structure”  of  turbulence,  that  is  it  will 
depend  for  instance  on  the  roughness  of  pipe  walls;  it  will  be  inde¬ 
pendent,  however,  of  the  physicochemical  properties  of  the  gas 
mixture,  and  especially  of  the  transformation  velocity  under  laminar 
conditions  S|Qm  . 


Very  few  experimental  data  are  available  for  comparison  with  this 

n 

theory*  Damkohler  carried  out  observations  on  Bunsen  burners,  which  in¬ 
dicate  an  increase  of  propagation  velocity  with  increasing  turbulence, 
but  do  not  allow  a  conclusive  quantitative  evaluation.  The  diffuse  char¬ 
acter  of  the  flame  front  under  turbulent  conditions  is  clearly  visible  in 

it 

Damkohler* s  flame  photographs.  Shchelkin  cites  work  on  combustion  in 
motors^  *  ,  as  experimental  evidence  for  his  theory;  the  increase  of 

propagation  velocity  with  the  number  of  revolutions  and  the  small  in¬ 
fluence  of  the  chemical  nature  of  the  fuel  observed  in  these  experiments 
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confirm  in  a  qualitative  way  the  theoretical  conclusions*  Obviously,  no 
quantitative  results  can  be  derived  because  of  the  lack  of  any  informa¬ 
tion  about  turbulent  conditions  in  engine  combustion.  Recent  observe- 
tions  of  Unger5  on  flame  propagation  in  tubes  show  the  increase  of 
propagation  velocity  with  increase  of  turbulence.  These  measurements 
were  carried  out  in  steel  tubes  in  vertical  position,  the  flame  propaga¬ 
ting  downward  against  the  gas  stream.  Figure  15  shows  a  plot  of  propa¬ 
gation  velocity  (calculated  as  the  sum  of  observed  flame  velocity  and 
average  gas  stream  velocity)  against  Reynolds  number  for  3  tube  diameters. 
Since  V  is  approximately  proportional  to  the  Reynolds  number,  the 
curves  do  not  agree  too  well  with  Shchelkin' s  hyperbolic  law,  but  they  show 
the  asymptotic  linear  relation  for  larger  values  of 


Shchelkin’ s  theory  was  criticized  and  a  modified  theory  was  proposed 
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recently  by  Wohl  .  The  expressions  proposed  by  this  author  ares 
(a)  for  JL  <  d  :  Sjurt,=  Sj0fn^/l  +■ 

with  the  approximations 

Sfurb  "  S|om+  k  V  for  small  V* 

Sturb5 


and 


Mom 


for  large  v* 


(  k  is  a  constant) 

(b)  for  Jl  >  d:  ,  S|urt' 

Slom 

with  S|gft,  "  \/ k  V*  S|om 

(c)  for  Jf>  d  ,  V  >  >  S|om 


Sturb*  ^  V 
Sturb 

as  approximate  solution  for  V  > 

^k>m 

_,*/3  «/3 

Sturb  "  k  V  S|om 
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In  concluding  this  chapter,  it  seems  hardly  necessary  to  point  out 
that  the  tydrodynamical  aspects  of  flame  propagation,  and  especially  the 
influence  of  turbulence,  are  highly  important  for  the  applications  of 
combustion ,  More  experimental  work  under  well-defined  conditions  of  tur¬ 
bulence  is  needed  to  evaluate  the  existing  theories.  From  the  theoreti¬ 
cal  point  of  view,  one  possible  effect  of  turbulence  on  flame  propaga- 

N 

tion  which  had  not  been  taken  into  account  by  Damkohler  and  by 
Shchelkin  is  a  direct  influence  on  the  course  of  the  chemical  reaction; 
for  instance  the  concentrations  in  the  reaction  zone  could  be  changed 
by  turbulence  in  such  a  way  that  the  average  life  time  of  active  par¬ 
ticles  would  be  modified.  Only  an  experimental  approach  to  this  problem 
seems  possible,  however,  because  of  tie  limited  knowledge  of  the  reaction 
mechanisms  involved. 


I «,  Theory  of  flames  formed  in  the  mixing  ?.one  of  fuel  and  air  (diffusion 
flames) 

In  the  preceding  chapters  it  had  been  assumed  tacitly  that  a  homo¬ 
geneous  combustible  mixture  of  fuel  and  air  (or  oxygen)  has  been  formed 
before  flame  propagation  takes  piece.  An  entirely  different  situation 
exists  when  the  mixing  of  fuel  gas  and  air  (or  oxygen)  and  combustion  take 
place  at  the  same  time  and  a  stationary  flame  is  formed  in  the  mixing  zone. 
In  the  important  technical  applications  of  this  type  of  combustion,  this 
mixing  is  always  a  result  of  turbulent  convection  and  diffusion.  However, 
only  the  case  of  true  diffusional  mixing  has  so  far  received  theoretical 


-  54  - 


treatment,  and  the  term  "diffusion  flame”  is  customary  in  this  case  . 
Since  diffusionel  mixing  occurs  much  more  siowly  than  the  chemical  re¬ 
action  could  proceed  at  the  usual  temperatures  and  pressures,  diffusion 
alone  is  the  rate-determining  factor  and  is  possible  to  treat  this  case 
without  introducing  any  assumptions  about  reaction  rates  and  mechanisms. 
On  the  other  hand  the  treatment  is  limited  to  a  determination  of  the 
shape  of  the  flame  and  no  information  about  the  rate  or  mechanism  of  re¬ 
action  can  be  expected.  This  problem  has  been  treated  by  Burke  and 
Schumann-^  (cf,  ref,  2,  pp  220-234)  for  the  cases  of  a  cylindrical  and 
of  a  flat  flame. 

Only  the  cylindrical  case  will  be  discussed  here.  Figure  16  shows 
the  experimental  arrangement  on  wnich  the  treatment  is  based.  The 
following  simplifying  assumptions  were  introduced  by  Burke  and  Schumann: 

(a)  The  flame  zone  coincides  with  the  locus  of  the  points  where  inter¬ 
diffusion  has  produced  a  stoichiometric  mixture. 

(b)  fuel  gas  and  air  flow  with  equal  velocity  V  ,  constant  over  the 
cross  section  as  well  as  along  the  flow. 

(c)  The  coefficient  of  ciffusion  is  constant, 

(d)  Diffusion  occurs  only  radially. 

Fortunately  the  effects  of  deviations  from  assumptions  (b)  and  (c) 
tend  to  compensate  each  other,  as  the  increase  of  temperature  in  the  flame 

*  It  is  hardly  necessary  to  point  out  that  this  case  is  entirely  different 
from  the  "di. ffusional  propagation"  discussed  on  page  34. 
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brings  about  at  the  same  time  an  increase  of  velocity  along  the  flow, 

increasing  the  length  of  the  flame,  and  on  increase  of  the  coefficient 

of  diffusion,  decreasing  the  flame  length.  In  order  to  simplify  the 

boundary  conditions  and  treat  the  problem  as  the  diffusion  of  a  single 

gas,  the  initial  concentration  of  oxygen  C2  is  introduced  as  a  negative 

C. 

concentration  of  fuel  — p  ,  where  f  is  the  number  of  oxygen  molecules 
which  combine  with  one  molecule  of  fuel.  By  setting  Y  :  V.t  the  diffusion 
equation  for  the  fuel  concentration  C  is  obtained  in  Uie  form 


Denoting  the  initial  fuel  concentrations  by  Cj  the  solution  must  ful¬ 
fill  the  following  conditions  (Figure  16) , 

c  =  for  Y  =  0,  0  =  r <  L 

c  =  _£t.  for  y  -  0  ,  L  r  R 

■  =0  for  r  50  and  r  -  R 

O  f 

Reference  is  made  to  the  original  paper  and  to  the  text  of  Lewis  and  von 
Elbe  for  the  details  of  the  mathematical  treatment.  By  setting  C  =  0, 
the  solution  gives  the  flame  shape  as  a  relation  bet?/een  r  and  Y  .  In 
the  case  of  an  underventilated  flame,  the  flame  front  will  curve  outward 
(Figure  16)  and  the  height  is  given  by  the  value  of  Y  for  r  =  R  ;  for  an 
overventilated  flame,  the  front  terminates  at  the  axis,  and  the  height  is 
the  value  of  Y  for  X  -  0  . 
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FIG. 

DIFFUSION  f 
(AFTER  BURKE  AND  SCI 


A  number  of  the  conclusions  reached  by  Burke  and  Schumann  are  de¬ 
rived  by  Jost  (ref  9.  pp.  207-216)  through  a  much  simpler  treatments 
He  points  out  that  the  assumptions  introduced  by  Burke  ana  Schumann  pre¬ 
clude  a  high  degree  of  accuracy,  so  that  his  treatment,  based  on  the 
same  assumptions,  should  be  adequate  enough.  The  depth  of  penetration  of 
the  air  into  the  fuel  gas  is  approximated  by  Jost  as  being  equal  to  the 
mean  square  displacement; 

X2  '  2  01  (66) 

2  9  Y 

By  identifying  X  with  r  and  t  with  there  results 


The  following  conclusions,  agreeing  with  those  of  Burke  and  Schumann, 
can  be  drawn  from  this  result •. 

(a)  Since  the  total  gas  flow  is  proportional  to  V  ffc  ,  the  flame  height 
Y  will  remain  unaltered  by  changing  the  dimensions  of  the  burner  (r) 

without  changir^  the  total  flow0 

(b)  The  height  Y  is  proportional  to  the  flow  velocity  V  and  inversely 
proportional  to  the  coefficient  of  dil fusion  D. 

(c)  For  constant  mass  flow,  a  change  in  pressure  will  not  affect  the 
height,  au  V  and  D  are  then  both  inversely  proportional  to  pressure. 

(d)  Substitution  of  an  inert  gas  for  part  of  the  combustible  gas  will 
shorten  an  overventilated  flame  and  lengthen  an  underventilated  one. 
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Jost’s  simple  treatment  allows  even  a  quantitative  estimation  of  the 
flame  height,  which  agrees  fairly  well  with  the  values  obtained  by 

Burke  and  Schumann  theoretically  and  experimentally.  As  en  example, 

* 

with  L  =  l/2  inch,  R  :  1  inch,  D  =  0.0763sq.  inch/sec.  and  V-0.610 
incty'sec.  the  rigorous  treatment  gives  Y  =  1.31  inch,  while  Jost's 
formula  leads  to 


0.61*0.5* 
"  2  x  0  0763 


1  inch 


Burke  and  Schumann  show  in  their  paper  the  good  agreement  with  experi¬ 
mental  results  achieved  by  their  theory „  A  still  better  agreement  was 
reached  hy  them  by  sampling  the  gas  composition  at  various  points  along 
the  axis  and  selecting  a  value  of  D  in  the  calculation  which  repro¬ 
duced  the  experimental  distribution  of  concentrations  as  closely  as 
possible. 


The  excellent  agreement  obtained  by  the  diffusion  theory  is  due,  as 
mentioned  before,  to  t?he  fact  that  the  rate  of  chemical  reaction  is 
large  compared  with  the  speed  of  the  oil  fusion  process.  On  the  other 
hand,  if  the  reaction  rete  is  small  compared  with  the  rate  of  diffusion, 
a  diffusion  flame  can  serve  as  a  meuns  to  investigate  reaction  rates 
(ref.  21,  p.  283,  ref.  40,  p.  15) 

Recently  Damkohler  and  Sander^  proposed  a  method  to  investigate 
fast  gas  reactions  under  approximately  isothermal  conditions,  based  on 
the  theory  of  diffusion  flames  of  Burke  and  Schumann.  The  arrangement  is 
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similar  to  the  one  shown  in  Figure  16  except  that  provision  is  made  for 
heating  the  gus  streams  before  entering  the  reaction  chamber  and  for  cool¬ 
ing  at  e  short  distance  downstream,  to  stop  the  reaction  and  allow  chemi¬ 
cal  anatysis  of  the  reaction  products  at  intermediate  stages.  The  relations 
derived  by  Burke  and  Schumann  were  found  to  hold  well  in  model  exper'znents 
with  02  and  C0g  as  interdif fusing  gases.  The  authors  point  out  that 
it  snculd  be  possible,  by  merely  changing  the  total  pressure,  and  keeping 
the  mass  flows  constant,  to  vury  the  reaction  time  witnout  changing  the 
temperature  and  concentration  distributions  in  the  reaction  chamber,  pro¬ 
vided  the  heat  of  reaction  can  be  neglected  compared  with  the  heat  content 
of  the  gases  The  length  of  the  reaction  chamber  is  chosen  so  small  that 
diffusion  to  the  wall  of  the  gas  entering  through  the  center  tube  is 
negligible,  precluding  the  possibility  of  a  catalytic  reaction  on  the  wall 

II 

surface.  No  actual  reaction  experiments  are  reported  by  Damkohler  and 
Sander,  so  that  the  value  of  this  application  of  diffusion  flames  for  the 
field  of  combustion  research  remains  to  be  proven. 

Even  with  combustion  in  a  homogeneous  mixture  diffusion  processes  of 
the  kind  discussed  in  this  chapter  may  under  certain  circumstances  play 
an  important  role.  Goldmanrt^*  suggested  that  the  fact  that  very  lean 
iydrogen-air  or  hydrogen-oxygen  mixtures  can  maintain  an  upward  propaga¬ 
ting  flame  could  be  explained  by  assuming  that  combustion  takes  place  in 
one  or  several  rising  gus  balls  which  receive  a  sufficient  supply  of  hy¬ 
drogen  by  diffusion.  This  accounts  also  for  the  fact  that  combustion  is 
less  complete  for  leaner  H2  mixtures.  A  modification  of  this  theory  has 
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been  proposed  by  Clusius,  Kolsch  and  Waldmann^.  In  their  combustion  ex¬ 
periments  with  lean  H?  —  Df  —  0Z  mixtures  these  authors  obtained  an  excess 
of  H  in  the  products  of  combustion  and  an  excess  of  02  in  the  remain¬ 
ing  mixture,  in  accordance  with  the  diffusion  theory  of  Goldmann,  Since 
they  could  not  observe  any  macroscopical  inhomogeneity  of  the  flame  front, 
they  suggested  that  the  combustion  zone  consists  of  small  filaments  per¬ 
pendicular  to  the  flame  front.  Combustion  takes  place  only  in  these 
filaments,  which  occupy  but  a  fraction  of  the  flame  surface,  and  hydrogen 
diffuses  into  the  filaments  and  maintains  thus  a  combustible  mixture. 

While  these  authors  did  not  actually  observe  this  structure,  some  observa¬ 
tions  of  a  filamentous  appearance  of  burner  flames  of  lean  hydrogen 
mixtures  and  of  rich  mixtures  of  higher  hydrocarbons  have  been  reported^. 

A  simple  explanation  of  the  way  in  which  the  flame  front  disinte- 
grates  into  these  filaments  has  been  given  by  Zeldovich  .  In  mixtures 
not  too  far  off  from  stoichiometric  composition  flame  propagation  is 
controlled  by  heat  conduction,  This  tends  to  stabilize  the  flame  front, 
because  the  regions  which  are  convex  toward  the  unburnt  gas  will  propa¬ 
gate  with  a  smaller  velocity  than  those  which  are  concave  (ref.  9,  p.125). 
In  very  lean  or  very  rich  mixtures,  however,  diffusion  will  be  the  con¬ 
trolling  process.  Since  diffusion  takes  place  in  the  direction  opposite 
to  the  propagation,  the  effect  will  be  contrary  to  that  of  heat  conduction: 
the  velocity  will  increase  in  the  convex  and  decrease  in  the  concave  parts, 
leading  ultimately  to  a  disruption  of  the  flame  front  (Figure  17). 
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r.  Conclusion  end  recommendations 

The  theory  of  flame  propagation  can  bo  subdivided  into  two  groups 
of  problems: 

(a)  study  of  the  transformation  velocity,  and 

(b)  hydrodynamics],  aspects  of  flame  propagation . 

The  theory  of  the  transformation  velocity  takes  into  account  the 
simultaneous  actions  of  heat  conduction,  diffusion  and  chemical  kinetics 
in  the  flame,  disregarding  any  effects  of  flew.  At  the  present  stage  of 
limited  knowledge  of  reaction  kinetics  the  treatments  are  based  on  hypo¬ 
thetical  assumptions  about  the  reaction  mechanisms.  So  far  the  theory  has 
been  developed  mainly  for  reactions  which  can  be  approximately  represented 

-  ]|y 

by  a  classical  kinetic  law,  W  =  S'6  j  even  in  this  simplest  case 

theoretical  values  do  not  agree  too  well  with  experimental  results.  Un¬ 
fortunately,  a  major  group  of  combustion  reactions  belongs  to  the  auto- 
catalytic  chain  type,  which  does  not  allow  such  an  approximation.  In 
addition  to  this  lack  of  definite  knowledge  about  reaction  kinetics,  which 
constitutes  the  main  obstacle  for  the  success  of  the  theory,  considerable 
difficulties  may  also  arise  for  the  treatment  cue  to  the  large  deviations 
from  thermal  equilibrium  likely  to  occur  in  the  combustion  zone.  It  is 
to  be  expected  that  new  methods  in  research  on  reaction  kinetics^  will 
provide  valuable  information  for  the  future  development  of  the  theory  of 
the  transformation  velocity. 
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The  hydrodynamical  treatment  of  flame  propagation  ha  a  so  far  given  a 
satisfactory  result  only  in  dealing  with  detonation,  and  with  the  simplest 
cases  of  slow  combustion,  namely,  the  spherically  symmetrical  flow  under 
conditions  of  constant  pressure  (soap  bubble)  and  constant  volume  (spherical 
bomb)  (cf.  ref„  1).  More  complicated  flow  conditions,  particularly  the 
important  cases  of  Bunsen  burner  flames  and  of  flame  propagation  in  pipes, 
have  been  treated  up  to  now  only  by  empirical  or  semi-cmpirical  methods* 

The  hydrodynamical  theory  of  one-dimensional  flow,  which  predicts  an 
upper  limit  for  the  flame  velocity  in  the  case  of  slow  combustion  (deflagra¬ 
tion),  has  not  been  able  to  account  for  the  fact  that  the  whole  range  from 
the  transformation  velocity  up  to  the  velocity  of  detonation  can  be  ob¬ 
served  under  suitable  experimental  conditions.  Several  different  flow 
effects  can  concur  in  producing  an  increase  of  the  propagation  velocity: 

(a)  turbulence  due  to  the  flow  in  the  unburnt  gas* 

(b)  turbulence  due  to  the  instability  of  the  flame  front  postulated 
by  Landau^,  and 

(c)  pressure  and  temperature  rise  due  to  the  compression  wave  gener¬ 
ated  in  the  unburnt  gas  by  the  combustion  process. 

The  theoretical  explanations  of  the  influence  of  hydrodynamical  pa¬ 
rameters  on  flame  propagation  usually  rest  on  the  basic  hypothesis  that 
the  transformation  velocity  is  not  affected  by  the  flow  conditions.  Ex¬ 
ceptions  to  this  assumption  are  made  only  for  the  case  of  turbulence 
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FIG.  18 

STROBOSCOPIC  PHOTOGRAPH  OF  A  BUNSEN  BURNER  FLAME 
DISTURBED  BY  AN  ALTERNATING  ELECTRICAL  FIELD. 

INNER  DIAMETER  OF  NOZ  2  LE  1.02  CM.  AVERAGE  FLOW 
VELOCITY  376  CM/SEC  425%  BY  VOLUME  NATURAL  GAS 
(NATOX  )  IN  AIR. 

2500  VOLT,  500  c.p.s.  APPLIED  BETWEEN  NOZZLE  AND 
SPHERE. 


t 


with  a  mixing  length  small  compared  to  the  thickness  of  the  flame  front, 
and  for  cases  in  which  the  velocity  of  the  unburnt  gas  relative  to  the 
flame  front  is  not  small  compared  with  the  velocity  of  sound.  In  all 
other  cases  the  influence  of  the  flow  parameters  is  assumed  to  be 
limited  to  an  increase  of  the  surface  of  the  reaction  zone. 

An  experimental  study  of  the  details  of  the  mechanism  by  which,  the 
surface  of  the  flame  is  increased  wr uld  furnish  valuable  information  for 
a  better  understanding  of  the  phenomena  encountered  in  flames  in  a  tur¬ 
bulent  medium.  It  appears  quite  hopeless,  however,  to  attempt  such  a 
study  under  conditions  of  actual  turbulence.  Preliminary  experiments 
carried  out  in  Cornell  Aeronautical  Laboratory^  indicate  that  one 
approach  to  this  problem  might  be  the  stroboscopic  observation  of  a 
flame  front  subjected  to  a  periodic  disturbance.  As  an  example,  Figure 
13  shows  the  stroboscopic  pattern  of  e  Bunsen  burner  flame  obteined  under 
the  influence  of  an  electrical  field  of  500  c  p.s.,  applied  between  the 
nozzle  and  the  spherical  electrode  visible  to  the  left  of  the  flame.  This 
method  should  provide  among  other  data  an  indication  of  the  upper  limit  of 
propagation  velocity  attainable  by  turbulence. 
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flow  velocity 
mass  flow 
density 

initial  density 
transformation  velocity 
heat  content 
vector  of  heat  flow 
reaction  rate 

heat  of  reaction  per  unit  mass  of  combustible 

coefficient  of  heat  conduction 

absolute  temperature 

specific  heat  at  constant  pressure 

initial  temperature 

ignition  temperature 

final  temperature,  temperature  of  combustion 
concentration  of  combustible 
initial  concentration  of  combustible 
pressure 

initial  pressure 

heat  of  reaction  per  unit  mass  of  total  gas 
mixture 

thickness  of  reaction  zone 
activation  energy 
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Re 
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gas  constant 

constants  1  n  express’ nr.  for  the  reaction  ra 

coefficient  of  diffusion 

coefficient  of  temperature  conduction 

kinematic  viscosity 

collision  coefficient 

mean  free  path 

mean  molecular  velocity 

tiiickness  of  heating  zone 

reel  part  of  (  ) 

characteristic  time,  reaction  Ume 

thicknes.*  of  fleme  zone 

turbulent  mixing  length 

instantaneous  value  of  fluctuation  velocity 
time  average  of  fluctuation  velocity 
coefficient  of  turbulent  exchange 
Reynolds  number 

flame  propagation  velocity  under  turbulent 
conditions 

fleme  propagation  velocity  under  laminar 
conditions 
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